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1.  The  technical  report  transmitted  herewith  represents  the  result  of 
Work  Unit  4A05  regarding  the  development  and  testing  of  a predictive 
model  for  marsh  productivity  and  succession.  This  work  unit  was  con- 
ducted as  part  of  Task  4A  (Marsh  Development)  of  the  Corps  of  Engineers' 
Dredged  Material  Research  Program  (DMRP) . Task  4A  is  a part  of  the 
Habitat  Development  Project  of  the  DMRP  and  is  concerned  with  develop- 
ing, testing,  and  evaluating  the  environmental,  economic,  and  engi- 
neering feasibility  of  using  dredged  material  as  a substrate  for  marsh 
development. 

2.  The  purpose  of  this  work  unit  was  to  develop  a simulation  model  of 
plant  growth  and  succession  on  salt  marshes.  The  model  was  created  by 
identification  of  the  biotic  and  physical  factors  that  control  plant 
growth  and  the  measurement  of  these  factors  under  field  conditions.  The 
model  appears  to  have  great  potential  for  use  in  establishing  criteria 
and  specifications  for  salt  marsh  development;  however,  additional 
parameterization  and  verification  will  be  required  before  the  system 
becomes  fully  operational. 

3.  Work  Unit  4A05  is  one  of  several  research  efforts  designed  by  the 
DMRP  to  accurately  document  marsh  productivity  and  the  factors  which 
influence  that  productivity.  Closely  related  work  units  are  4A04A1  and 
4A04B,  which  address  the  productivity  of  minor  marsh  species  along  the 
Atlantic  and  Gulf  coasts,  respectively;  4A04A2,  which  deals  with  marsh 
plant  substrate  selectivity  and  underground  biomass  production;  and 
4A20,  a less  intensive  effort  that  will  provide  a general  evaluation  of 
salt  marsh  productivity  on  the  Pacific  coast  of  the  United  States. 
Additional  supportive  and  comparative  data  will  be  forthcoming  with  the 
final  analysis  of  the  results  of  field  studies  at  Windmill  Point, 
Virginia,  (4A11);  Buttermilk  Sound,  Georgia,  (4A12) ; Apalachicola, 
Florida,  (4A19) ; Bolivar  Peninsula,  Texas,  (4A13) ; Pond  No.  3,  Cali- 
fornia, (4A18);  and  Miller  Sands,  Oregon,  (4B05) . Together  these 
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research  products  provide  the  Corps  with  a comprehensive  basis  for  sound 
management  decisions  regarding  habitat  development  on  dredged  material 
and  disposal  in  natural  marsh  habitats. 


//JOHN  L.  CANNON 

Colonel,  Corps  of  Engineers 
Commander  and  Director 
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spicata  were  more  variable,  but  where  regular  flooding  occurred  there  was  a 
strong  positive  correlation  with  !>.  patens  abundance  and  negative  correlation 
with  D.  spicata. 

Records  from  the  continuous  monitoring  study  reinforce  the  concept  of 
tidal  Inundation  acting  as  a master  variable.  Power  spectrum  analysis  of 
continuous  research  of  temperature  and  pH  in  the  sediments  revealed  that  in 
the  S_.  alterniflora  zone  the  sediment  surface  temperature  showed  the  major 
power  peak  to  be  0.089  cycle/hour,  which  is  the  tidal  frequency  for  this  area. 
In  the  higher  £.  patens  - D.  spicata  zone  the  frequency  with  the  largest  power 
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radiation.  The  analysis  of  pH  data  showed  a large  influence  in  the  low  £. 
alterniflora  at  the  tidal  frequency  but  not  in  the  high  marsh. 

y*  A model  was  developed  which  depicted  plant  growth  and  succession  across 
a transect.  The  model  was  driven  by  solar  radiation,  temperature,  soil 
salinity,  and  tidal  inundation.  The  model  showed  plant  growth  controlled 
primarily  by  tidal  inundation,  with  Influences  felt  from  radiation  and  tem- 
perature, but  little  influence  due  to  salinity.  The  model  appears  to  have 
great  utility  for  future  use  but  presently  is  limited  in  application  by 
incomplete  parameterizations.  a 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAOE(Vh»n  Dmtm  Bnfnd) 


I 


EXECUTIVE  SUMMARY 

Rapid  disappearance  and  deterioration  of  wetlands  as  a conse- 
quence of  increasing  human  utilization  coupled  with  growing  recogni- 
tion by  the  scientific  community  of  the  ecological  importance  of  these 
areas  Indicate  the  need  to  Investigate  ecologically  sound  alternatives 
for  using  dredged  material.  Attempts  to  recolonize  disposal  areas 
after  deposition  have  thus  far  met  with  mixed  success.  In  attempting 
to  establish  a desired  community,  two  options  are  available:  (1)  create 
environmental  conditions  favorable  for  the  colonization  of  the  desired 
species,  or  (2)  introduce  a species  that  is  compatible  with  the  present 
environmental  conditions  and  allow  natural  succession  to  produce  the 
desired  species  or  community.  In  order  to  achieve  the  desired  ecosys- 
tem by  either  of  the  two  above  approaches,  the  dominant  controlling 
forces  on  the  communities  must  be  known  and  described. 

Three  parallel  studies  were  conducted  with  the  development  of  a 
simulation  model  of  plant  growth  and  succession  on  a salt  marsh  being 
the  ultimate  objective.  The  first  component  was  concerned  with  deter- 
mining the  biotic  and  chemical  parameters  present  at  a series  of 
marshes.  The  second  component  was  concerned  with  obtaining  continuous 
records  of  several  physical  and  chemical  parameters  within  several 
vegetation  zones.  The  third  component  of  the  project  was  the  development 
of  the  simulation  model. 

Field  ;tudies  showed  that  Spartina  altemiflora  exhibited  a 
strong  positive  correlation  with  elevation  and  tidal  inundations  at  the 
two  major  study  sites,  Taskinas  Creek  and  Bennett's  Creek,  in  Virginia. 
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S.  alterniflora  showed  a negative  correlation  with  soil  salinity  and 
significant  positive  correlation  with  available  iron.  Responses  of 
j>.  patens  and  Distichlis  spicata  were  more  variable,  but  where  regular 
flooding  occurred  there  was  a strong  positive  correlation  with  S. 
patens  abundance  and  negative  correlation  with  D.  spicata. 

Records  from  the  continuous  monitoring  study  reinforce  the  con- 
cept of  tidal  inundation  acting  as  a master  variable.  Power  spectrum 
analysis  of  continuous  research  of  temperature  and  pH  in  the  sediments 
revealed  that  in  the  S..  alterniflora  zone  the  sediment  surface  tempera- 
ture showed  the  major  power  peak  to  be  at  0.089  cycle/hour,  which  is 
the  tidal  frequency  for  this  area.  In  the  higher  S^.  patens  - D. 
spicata  zone  the  frequency  with  the  largest  power  component  is  that 
of  nearly  1 cycle/day  corresponding  to  the  input  of  solar  radiation. 

The  analysis  of  pH  data  showed  a large  influence  in  the  low  S^.  alterni- 
flora at  the  tidal  frequency  but  not  in  the  high  marsh. 

A model  was  developed  which  depicted  plant  growth  and  succession 
across  a transect.  The  model  was  driven  by  solar  radiation,  tempera- 
ture, soil  salinity,  and  tidal  inundation.  The  model  showed  plant 
growth  controlled  primarily  by  tidal  inundation,  with  influences  felt 
from  radiation  and  temperature,  but  little  influence  due  to  salinity. 
The  model  appears  to  have  great  utility  for  future  use  but  presently 
is  limited  in  application  by  Incomplete  parameterlzations. 
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PREFACE 


This  report  presents  the  results  of  an  investigation  of  factors 
influencing  plant  zonation  and  succession  on  salt  marshes.  The  study 
was  supported  by  the  U.  S.  Army  Engineer  Waterways  Experiment  Station 
(WES),  Habitat  Development  Project,  Dredged  Material  Research  Program, 
Vicksburg,  Mississippi,  under  contract  No.  DACW39-73-C-0121  to  the  Uni- 
versity of  Virginia,  Charlottesville,  Virginia. 

The  project  was  directed  by  Joseph  C.  Zieman  and  William  E.  Odum 
of  the  Department  of  Environmental  Sciences  of  the  University  of 
Virginia.  Field  Sampling  and  laboratory  analysis  were  performed  by 
Mark  L.  Gleason,  Thomas  G.  Wolaver,  Jeffrey  E.  Drifmeyer,  and  Michael  B. 
Robblee . 

The  study  was  conducted  under  Dr.  Kent  Thornton,  Contract  Manager, 
and  Dr.  Hanley  K.  Smith,  Habitat  Development  Project  Manager,  WES.  The 
study  was  under  the  general  supervision  of  Dr.  John  Harrison,  Chief, 
Environmental  Effects  Laboratory,  WES. 

The  Directors  of  WES  during  the  study  and  the  preparation  of  this 
report  were  COL  G.  H.  Hilt,  CE,  and  COL  J.  L.  Cannon,  CE.  Technical 
Director  was  Mr.  F.  R.  Brown. 
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Depicts  the  simulation  at  the  6-ft 
elevation  above  mlw  on  the  marsh. 

(1.83-m) 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


The  U.  S.  customary  units  of  measurement  used  in  this  report  can  be 
converted  to  metric  units  as  follows: 


Multiply 

By 

To  Obtain 

microns 

0.000001 

meters 

inches 

2.54 

centimeters 

feet 

0.3048 

meters 

acres 

4046.856 

square  meters 

pounds  (force)  per 
square  incfc 

6894.757 

pascals 

miles  (U.  S. 
statute) 

1.609344 

kilometers 
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PART  I : INTRODUCTION 

1.  In  recent  years  the  value  of  tidal  wetland  as  nursery  areas,  as 
buffers  against  shoreline  erosion,  and  as  sources  of  energy  for  shell- 
fish and  commercial  and  sport  fisheries  has  become  widely  recognized 
(Schelske  and  Odum,  1971;  Teal,  1962;  Odum  and  de  la  Cruz,  1967; 

Cooper,  1969;  Williams  and  Murdoch,  1969;  Keefe,  1972;  Odum,  Zieman 

and  Heald,  1972).  Estuaries  are  among  the  most  productive  natural  eco- 
systems on  earth,  and  serve  as  vital  recreational,  transportation, 
industrial,  commercial,  and  population  sites  (Allen,  1964). 

2.  With  accelerating  human  uses  of  estuaries  have  come  altera- 
tions and  modifications  of  the  natural  ecosystem.  Chief  among  these 
is  dredging  and  dredged  material  disposal.  Dredging  and  filling  has 
altered  over  7 percent  of  the  total  estuarine  habitat  in  the  United 
States;  570,000  out  of  8 million  acres  have  been  lost  (Marshall,  D.  E., 

1970).  Along  the  northeastern  Atlantic  coast,  dredging  is  the  predominant 
means  of  wetland  alteration  (Marshall,  D.  E.,  1970;  Schmidt,  1966).  For 
example,  of  the  45,000  acres  of  wetlands  between  Maine  and  Delaware 
destroyed  from  1955  to  1964,  fully  one-third  was  lost  to  dredged  mate- 
rial disposal  (Schmidt,  1966).  During  this  same  period.  South  Carolina 
lost  10  percent  of  its  total  wetland  acreage  to  dredged  material  dis- 
posal, while  Florida  filled  60,000  acres,  or  7.5  percent,  of  its  most 
important  estuarine  areas  (Teal  and  Teal,  1969;  Marshall,  D.  E.,  1970). 

Along  the  Gulf  Coast  dredging  and  fill  operations  have  also  impacted 

* A table  of  factors  for  converting  U.S.  customary  units  of  measure- 
ment to  metric  (SI)  units  can  be  found  on  page  12. 
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estuaries;  Texas  lost  55,000  acres  of  shallow  marine  bay  habitat  and 
23,000  acres  of  brackish  marsh  (Chapman,  1967). 


' 


3.  The  biological  effects  of  dredging  and  disposal  are  often 
deleterious,  as  has  been  discussed  by  Thompson,  1961;  Chapman,  1967; 
and  Odum,  1970.  Excessive  sedimentation,  increased  turbidity  and 
alteration  of  circulation  patterns  may  adversely  affect  benthic  popula- 
tions, photosynthetic  activity,  and  energy  flow.  For  example  in  Boca 
Ciega  Bay,  Florida,  3,500  acres  of  dredged  material  deposition  severely 
reduced  annual  production  of  seagrasses  (by  25,851  metric  tons), 
infauna  (by  1,091  metric  tons),  and  fisheries  (by  73  metric  tons) 
(Marshall,  D.  E.,  1970).  While  not  an  exhaustive  list,  the  above  show 
that  dredging  operations  have  the  potential  to  alter  significant  portions 
of  existing  wetlands. 

4.  Rapid  disappearance  and  deterioration  of  wetlands  as  a con- 
sequence of  burgeoning  human  utilization,  coupled  with  growing  recogni- 
tion by  the  scientific  community  of  the  ecological  importance  of  these 
areas,  indicate  the  need  to  investigate  ecologically  sound  alterna- 
tives for  using  dredged  material.  In  particular,  the  possibility  of 
using  dredged  material  to  build  new  marshes  has  been  suggested 
repeatedly  (Smith,  1942;  Larimer,  1968;  Teal  and  Teal,  1969;  Woodhouse 
et  al.,  1972;  Connecticut  State  Board  of  Fisheries  and  Game,  undated). 
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PART  II:  PROJECT  CONCEPT  AND  DEVELOPMENT 


5.  Understanding  of  marsh  plant  physiology  and  ecology  is  pre- 
sently insufficient  to  offer  design  criteria  for  dredged  material  dis- 
posal that  would  facilitate  rapid  recolonization  of  marshes.  Attempts 
to  recolonize  disposal  areas  after  deposition  have  met  with  mixed  suc- 
cess (Windom,  1972;  Woodhouse  et  al.,  1972;  Vittor,1971;  Larimer,  1968). 
Ecological  succession  is  an  orderly  process  of  community  development 
that  involves  changes  in  species  composition  and  community  processes 
over  time.  In  marsh  areas,  the  sequence  of  plant  species  is  well  docu- 
mented, but  the  processes  and  driving  forces  behind  these  changes  are 
poorly  understood.  As  one  community  (an  assemblage  of  populations 
within  a given  area)  grows,  it  modifies  its  environment,  creating  con- 
ditions favorable  for  a subsequent  community.  The  environment  of  a 
mature  community  may  resemble  that  of  its  early  successional  predeces- 
sors only  slightly — a fact  that  is  often  overlooked  by  those  who  attempt 
to  establish  particular  communities  artificially. 

6.  In  attempting  to  establish  a desired  community,  two  options 
are  available:  (1)  create  the  environmental  conditions  favorable  for 
the  colonization  of  a desired  species,  or  (2)  introduce  a species  that 
is  compatible  with  the  present  environmental  conditions  and  allow 
natural  succession  to  produce  the  desired  species  or  community.  In 
order  to  achieve  the  desired  ecosystem  by  either  of  the  above 
approaches,  the  dominant  controlling  forces  on  the  communities  must  be 
known  and  described.  In  addition,  the  second  procedure  requires  know- 


ledge of  the  direction  in  which  succession  will  proceed. 
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7.  A comprehensive  understanding  of  the  development  of  estuarine 


marshes  requires  knowledge  of  the  physical  and  chemical  parameters 
affecting  marsh  plants.  The  ways  in  which  sediment  characteristics 
such  as  organic  content,  sediment  texture,  and  mineralogy  condition 
the  habitat  of  marsh  plants  have  been  elucidated  by  previous  workers 
(Richards,  1954;  Fleming  and  Alexander,  1961;  Windom,  1972). 

8.  According  to  the  classification  of  Sanders  (1969),  marsh  eco- 
systems are  controlled  primarily  by  physical  factors.  The  ecological 
stresses  imposed  by  the  broad  variability  of  physical  factors  reduce 
the  relative  influence  of  biotic  competition  as  a control  of  species 
distribution.  Since  relatively  few  plants  and  animals  can  withstand 
rapid  tidal,  salinity,  and  temperature  fluctuations,  species  diversity 
in  tidal  marshes  tends  to  be  low.  In  a study  of  saline  and  brackish 
marshes  in  Louisiana,  Chabreck  (1972)  found  that  in  the  low  marsh  zone, 
dominated  by  Spartina  altemiflora  and  strongly  influenced  by  tides, 
three  species  accounted  for  over  70  percent  of  the  higher  plants 
encountered . 

9.  The  tools  of  simulation  modeling  which  are  beginning  to  come 
into  use  in  environmental  studies  offer  great  potential  for  the  future 
(Patten,  1971,  1972).  However,  most  attempts  to  apply  modeling  tech- 
niques to  ecological  problems  have  suffered  from  several  crucial 
deficiencies.  First,  to  be  a useful  decision  tool,  a model  must  be 
designed  to  answer  specific  questions.  Thus  far,  few  ecological  models 
have  been  so  structured;  rather,  most  attempt  to  describe  the  "normal 
behavior"  of  a particular  system,  and  have  thus  been  largely 
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descriptive.  Second,  to  properly  develop  a simulation  model  of  an  eco- 
system it  is  necessary  to  understand  component  interactions  and  the  bio- 
logical theory  behind  these  interactions.  Third,  for  proper  implementa- 
tion the  best  conceived  model  still  requires  reliable  dynamic  informa- 
tion, such  as  flux  rates.  This  information  has  seldom  been  measured  in 
the  past  and  generally  is  not  available  in  the  literature.  For  best 
results,  the  input  parameters  should  be  measured  simultaneously  and  con- 
tinuously. 

10.  At  present  the  difficulty  of  collecting  accurate  dynamic 
field  data  and  ignorance  of  the  biological  functioning  of  the  marsh  eco- 
system are  more  serious  hindrances  to  ecological  modeling  than  misappli- 
cations of  systems  analysis  and  modeling  techniques.  One  of  the 
toughest  problems  faced  by  ecological  modelers  today  is  in  the  gather- 
ing of  accurate  environmental  data  at  closely  spaced  intervals  of  time. 
For  years,  physical  scientists  have  recorded  data  at  exceptionally  high 
acquisition  rates;  biologists,  however,  have  been  restricted  to  sam- 
pling field  sites  on  a weekly  or  monthly  basis,  distorting  interpreta- 
tion of  environmental  parameters  that  vary  widely  and  rapidly. 

11.  Simulation  modeling  is  essentially  a cyclical  process  cen- 
tered around  the  design  and  conceptualization  of  the  model,  the  formu- 
lation process,  the  parameterization  of  this  model,  and  field  studies 
to  verify  the  model.  Following  the  completion  of  one  cycle,  the  informa- 
tion gained  is  used  to  modify  the  original  concept  and  resulting  model. 
The  process  can  be  repeated  until  satisfactory  results  are  produced. 

12.  The  ultimate  objective  of  this  project  was  to  produce  an 


- 
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ecological  simulation  model  capable  of  predicting  vegetation  succession 
on  marshes  and  dredged  material  islands.  Three  parallel  studies  were 
pursued.  The  first  component  was  concerned  with  determining  the  biotic 
and  chemical  parameters  in  a series  of  marshes.  The  objectives  were  to 
determine  the  range  of  selected  environmental  parameters  that  dominant 
species  of  marsh  plants  tolerate,  and  second  to  determine  the  seasonal 
variation  of  selected  physicochemical  parameters  and  the  response  of 
certain  biotic  compounds.  A series  of  marshes  from  New  Jersey  to 
Georgia  were  visited  and  surveyed,  but  intensive  effort  was  concentra- 
ted at  two  sites  in  Virginia.  The  results  of  this  study  are  discussed 
in  Part  V. 

13.  The  second  study  was  concerned  with  obtaining  continuous 
records  of  several  physical  and  chemical  parameters  within  several 
vegetative  zones.  The  objective  of  this  study  was  to  characterize  the 
magnitude  of  variation  of  edaphic  factors  in  the  major  marsh  communi- 
ties. This  involved  the  development  of  a continuous  recording  data- 
monitoring  system  described  in  Appendix  A and  the  data  and  results  are 
discussed  in  Part  VI. 

14.  The  third  study  of  the  project  was  the  development  of  a 
simulation  model  of  plant  growth  and  succession  on  the  marsh.  The 
iritial  models  were  developed  based  on  information  available  in  the 
literature  on  the  functioning  of  estuarine  salt  marshes.  As  the  pro- 
ject progressed,  information  was  added  from  field  studies  and  contact 
with  other  scientists  in  the  field.  The  latter  input  was  aided  by  sev- 
eral meetings  on  marsh  growth  and  reclamation  at  the  Waterways 


Experiment  Station,  Vicksburg,  Mississippi.  Part  VII  contains  the 
development  of  the  marsh  succession  model. 


PART  III:  STUDY  SITES 


General  Description 

15.  Field  research  was  conducted  Intensively  at  two  sites  in 
Virginia  as  described  below.  In  addition,  surveys  were  conducted  at  a 
variety  of  marshes  from  New  Jersey  to  Georgia.  The  data  from  the  marsh 
surveys  are  presented  in  Appendix  C. 

16.  Taskinas  Creek  Marsh  is  situated  adjacent  to  the  York  River 
(lat.  37°25',  long.  76°43'),  16  km  southeast  of  the  city  of  West  Point. 
The  marsh  and  much  of  its  watershed  are  included  within  York  River 
State  Park  (Figure  1 ) . The  second  marsh  site  is  located  on  the  east 
bank  of  Bennett's  Creek,  a tributary  of  the  Nansemond  River,  about 

300  m south  of  U.S.  Route  17  (lat.  36°51',  long.  76°28'30"),  8 km  west 
of  Portsmouth.  The  marsh  and  adjacent  uplands  are  privately  owned. 

17.  Although  they  are  approximately  60  km  apart,  Taskinas  Creek 
Marsh  and  Bennett's  Creek  Marsh  are  alike  in  several  respects.  Both 
are  small  (less  than  37  ha)  pocket  marshes  adjoining  the  tidal  tribu- 
taries of  Chesapeake  Bay.  Their  average  tidal  range  (approximately 

1 m)  differs  by  less  than  5 cm.  Salinity  of  the  creek  water  ranges 
from  zero  to  approximately  15  percent  at  both  sites.  The  distribution 
of  vegetational  communities  is  also  similar.  Spartina  altemif lora 
dominates  the  lower  margin  of  each  marsh  along  the  creek  banks  while  a 
Spartina  patens  - Distichlis  spicata  association  occupies  the  adjacent 
high  marsh  zone.  Scirpus  olney , and  Scirpus  robustus  predominate 
where  the  marsh  borders  the  uplands.  The  average  width  of  each  marsh 
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from  creek  bank  to  upland  Is  approximately  100  meters. 

Marsh  Development  History 

18.  Abundant  evidence  suggested  that  vegetation  distribution  can 
change  extensively  as  a marsh  develops  (Redfield  and  Rubin,  1962; 
Redfield,  1967;  Redfield,  1972;  Newman  and  Rusnak,  1965;  Bloom,  1964; 
KnighL,  1934).  The  frequently  stated  hypothesis  that  plant  species 
distribution  was  somehow  intimately  influenced  by  elevation  of  the 
marsh  surface  relative  to  tidal  changes  in  water  level  (Redfield,  1972; 
Johnson  and  York,  1915;  Adams,  1963)  implied  that  a change  in  elevation 
will  alter  the  species  distribution.  Consequently,  the  history  of  plant 
distribution  in  a tidal  marsh  reflected  a dynamic  balance  among  factors 
that  influence  elevation:  accretion  of  sediment  on  the  marsh  surface, 
compaction  of  sediments  after  burial,  eustatic  changes  in  sea  level,  and 
subsidence  of  uplifting  of  the  land  mass  (Bloom,  1964;  Newman  and  Rusnak, 
1965;  Redfield,  1972;  Redfield  and  Rubin,  1962).  Reconstructing  the  evo- 
lution of  a marsh  can  furnish  a useful  context  for  understanding  present 
patterns  of  plant  distribution;  accordingly  a brief  study  was  undertaken 
to  outline  the  developmental  histories  of  Taskinas  Creek  Marsh  and 
Bennett's  Creek  Marsh. 

19.  The  peaty  sediments  underlying  the  marsh  surface  retain  ves- 
tiges of  the  environment  that  prevailed  when  the  sediments  were  depos- 
ited (Redfield,  1967).  In  order  to  gain  a general  picture  of  environ- 
mental change  over  time,  peat  cores  were  sampled  from  high  and  low 
marsh  stations  at  each  marsh,  using  an  Oakfield  Apparatus  Company 
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sectioned  peat  corer,  to  a maximum  depth  of  8.5  m.  Sixteen-centimeter 


core  sections  were  retained  at  1.25-m  depth  intervals  for  the  following 
laboratory  analyses:  percent  moisture,  percent  organic  matter,  and 
plant  species  identification. 

20.  Core  sections  were  analyzed  for  moisture  and  organic  matter 
as  outlined  under  "Laboratory  Procedures"  in  Part  IV.  Subsamples  of 
each  section  were  washed  through  three  sediment  screens,  of  mesh 
sizes  2 mm,  250  pm,  and  63  pm,  respectively.  The  fraction  retained 

on  each  screen  was  examined  under  a dissecting  microscope;  when  possi- 
ble, root  and  stem  remains  were  identified  by  species.  Each  fraction 
was  also  described  qualitatively,  noting  abundance  of  organic  matter, 
presence  of  foraminiferan  tests,  type  of  inorganic  sediment,  and  other 
observable  features. 

21.  Approximate  stratigraphic  sequences  in  each  marsh  have  been 
reconstructed  from  the  data.  Figure  2 depicts  cross  sections  on  low 
marsh  to  high  marsh  transects  at  Taskinas  Creek  Marsh  and  Bennett's 
Creek  Marsh,  respectively.  Solid-line  borders  between  adjacent  strata 
connect  points  for  which  data  exist;  dashed  lines  represent  more  specu- 
lative extrapolations. 

22.  The  deepest  stratum  sampled  at  Taskinas  Marsh  (Figure  2) 
appeared  to  be  sedge  peat,  characterized  by  sparsely  scattered  sedge 
stems.  Overlying  the  sedge  peat  was  a relatively  thick  wedge  of  fine- 
textured  clay  silt,  nearly  devoid  of  identifiable  detritus  particles. 
Above  the  mud  layer  was  a thick  wedge  of  dense,  fibrous  high  marsh  peat 
composed  predominantly  of  j>.  patens  and  D.  spicata  rhizome  remains, 
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which  extended  uninterrupted  to  the  marsh  surface  in  the  present-day 
high  marsh.  The  most  recent  deposit  was  a relatively  shallow  (.7-1.0  m) 
lens  of  low  marsh  (S.  alternif lora)  peat  covering  the  modern  intertidal 
zone. 

23.  The  stratigraphic  sequence  at  Bennett's  Marsh  (Figure  2 ) 
was  similar  to  the  Tasklnas  profile,  but  with  several  important  distinc- 
tions. Terrestrial  sediments  found  at  a depth  of  6.8  m at  station  #4  in- 
dicated that  the  total  thickness  of  marsh  sediments  may  be  considerably 
less  than  at  Taskinas  Creek — at  least  1.9  m less.  Ac  at  Taskinas  Marsh,  a 
layer  of  sedge  remains  was  covered  by  a thick  deposit  of  silty  clay.  How- 
ever, low  marsh  deposits  originated  considerably  deeper  at  Bennett's  than 
at  Taskinas,  while  high  marsh  peat  at  Bennett's  that  underlaid  the  modern 
high  marsh  surface  was  comparatively  shallow  (less  than  0.6  m) . 

24.  Lacking  radiocarbon  analyses  of  the  peat  samples,  it  was 
difficult  to  accurately  establish  the  age  of  each  stratum.  Neverthe- 

1 

less,  a general  sequence  of  events  was  deduced  by  the  principle  of 
superposition  (Leet  and  Judson,  1965). 

25.  Scattered  studies  of  marshes  bordering  the  East  Coast  of  the 
United  States  suggested  that  contemporary  marshes  were  of  relatively 
recent  origin.  Bloom  (1964)  estimated  a maximum  age  of  approximately 


7000  years  before  present  (B.P.)  for  peat  at  the  basement  of  a 
Connecticut  salt  marsh.  Redfleld's  (1967,  1972)  Cape  Cod  marsh  dates 
back  to  4000  years  B.P.,  while  the  Wachapreague  Marshes  of  Virginia's 
eastern  shore  appeared  about  5000  years  ago  (Newman  and  Rusnak,  1965). 
It  seemed  reasonable,  then,  to  infer  that  Taskinas  Creek  Marsh  and 


Bennett's  Creek  Marsh  originated  within  the  last  4000  to  7000  years. 
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26.  Based  on  core  samples,  the  earliest  macrophyte  community  at 
both  marshes  was  apparently  a brackish-water  association  characterized 


by  sedges  such  as  Scirpus . Bloom  (1964)  found  a similar  association 
in  Connecticut  marsh  deposits.  These  early  sedge  marshes  may  have 
formed  as  narrow  fringes  separating  the  uplands  from  adjacent  shallow 
bays.  As  sea  level  rose  relative  to  the  land,  the  fringing  marshes 
continued  to  colonize  the  borders  of  the  drowning  uplands,  while  their 
seaward  edges  were  reclaimed  by  the  bays.  Even  today,  the  most  abun- 
dant stands  of  Scirpus  in  Virginia  salt  marshes  congregate  along  the 
upland  margins.  The  association  of  a sedge  zone  fringing  shallow  bays 
has  also  been  identified  in  Connecticut  (Bloom,  1964),  and  elsewhere  in 
Virginia  (Newman  and  Rusnak,  1965) . 

27.  The  histories  of  Taskinas  Marsh  and  Bennett's  Marsh  appeared 
to  diverge  with  the  initiation  of  a high  marsh  (S.  patens  and  D.  spi- 
cata)  community  at  Taskinas,  at  a maximum  depth  of  about  5.5  m below 


the  present  marsh  surface.  The  data  suggested  that  the  transition 
from  shallow  bay  mud  to  fibrous  high  marsh  peat  occurred  relatively  rap- 
idly, perhaps  even  without  an  intervening  low  marsh  stage  (Redfield, 
1972).  Expansion  of  the  high  marsh  may  have  coincided  with  the  slow- 
down in  rate  of  relative  sea  level  rise  that  occurred  along  the  East 
Coast  around  2000-3000  years  B.P.,  facilitating  uplift  of  the  marsh 
surface  by  sediment  accretion  (Redfield  and  Rubin,  1962;  Newman  and 
Rusnak,  1965;  Bloom,  1964).  The  high  marsh  community  spread  across 


Taskinas  Marsh,  until  the  relatively  recent  appearance  of  alterni- 
flora  along  the  creek  banks.  Organic  remains  in  the  upper  0.9  m of 
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sediment  indicated  progressive  expansion  of  this  low  marsh  species  rear- 
ward across  previous  high  marsh.  Furthermore,  Taskinas  high  marsh  was 
dotted  with  Iva  frutescens,  a shrub  generally  associated  with  a more 
elevated  habitat  (Kurz  and  Wagner,  1957).  Almost  all  of  the  I. 
frutescens  were  dead.  Coupled  with  the  invasion  of  J>.  alterniflora 
over  the  high  marsh,  these  observations  suggested  that  Taskinas  Marsh 
is  presently  undergoing  a drowning  phase;  that  is,  the  rise  in  relative 
sea  level  is  outstripping  the  vertical  accretion  of  sediment,  and  retro- 
gressive succession  (Kurz  and  Wagner,  1957)  appears  to  be  occurring. 

28.  The  mud  stratum  at  Bennett's  Marsh  extended  upward  consider- 
ably closer  to  the  surface  than  at  Taskinas  Marsh.  Although  relative 
ages  of  similar  strata  in  each  marsh  could  not  be  determined  without 
radio-carbon  techniques,  it  seemed  plausible  that  an  open-water  environ- 
ment persisted  until  more  recently  at  Bennett's  Marsh,  particularly  in 
the  zones  now  occupied  by  high  marsh  species. 

29.  It  appeared  that  a j>.  alterniflora  community  occupied  the 
intertidal  zone  at  Bennett's  long  before  a J5.  patens  - D.  spicata 
association  became  established,  since  the  high  marsh  peat  layer  is  much 
shallower  than  the  low  marsh  peat  (Figure  2) . The  presence  of  low 
marsh  and  high  marsh  associations  above  the  mud  stratum  suggested  chat 
Bennett's  Marsh  has  recently  been  undergoing  an  accreting  phase;  that 
is,  the  vertical  accumulation  of  sediment  is  surpassing  the  rise  in 
relative  sea  level.  However,  the  data  were  insufficient  to  indicate 
whether  high  marsh  is  now  succeeding  low  marsh,  or  vice  versa,  or 
whether  the  balance  is  stable. 
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30.  In  summary,  the  marsh  stratigraphy  investigation  afforded 


t 


the  following  perceptions: 

(a)  The  development  of  Taskinas  Marsh  and  Bennett's  Marsh 
paralleled  the  general  patterns  outlined  for  other  East  Coast  salt 
marshes  (Redfield,  1967,  1972;  Newman  and  Rusnak,  1965;  Bloom,  1964; 

Hill  and  Shearin,  1970). 

(b)  Taskinas  Marsh  is  presently  undergoing  a drowning  phase, 
as  sea  level  has  risen  relative  to  the  surface  of  the  marsh. 

(c)  Bennett's  Marsh  has  recently  been  experiencing  an 
accreting  episode;  i.e.,  the  marsh  surface  has  been  elevated  relative 
to  sea  level. 
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PART  IV : METHODS 

Field  Procedures 

31.  After  the  results  of  preliminary  studies  were  evaluated,  a 
final  transect  design  was  adopted  at  the  two  intensive  study  sites, 
Taskinas  Creek  and  Bennett's  Creek.  Three  parallel  belt  transects  were 
established  in  each  marsh,  extending  across  the  ecotone  from  low  marsh 
to  high  marsh.  The  transects  (Figure  3)  were  8 m longer  at  Bennett's 
Marsh  than  at  Taskinas  Marsh,  and  thereby  encompassed  two  more  sample 
zones  in  order  to  accommodate  Bennett's  broader  distribution  of  vegeta- 
tive zones. 

32.  Within  each  transect  were  equidistant  "sample  zones",  4-m  x 
1-m  areas  from  which  all  plant  and  sediment  samples  were  taken  (Figure 
3) . Samples  were  gathered  randomly  within  each  zone  by  means  of  ran- 
domized assignment  of  numbers  of  the  four  adjacent  1-m2  quadrats  in 
each  sampling  zone. 

33.  Plant  samples  and  sediment  cores  were  taken  at  four  times 

during  the  spring  and  summer  of  1975:  5-7  April,  29-30  May,  21  July, 

and  13  September.  On  each  occasion,  a single  plant  sample  (shoots 
only)  and  three  sediment  cores  were  taken  from  each  sample  zone.  Qua- 
drats harvested  were  0.1  m2  in  zones  where  S^.  alternif lora  appeared  to 
be  prevalent,  and  0.04  m2  where  S..  patens  and  D.  spicata  predominated. 
Cores  were  taken  with  a brass  coring  apparatus  (6.35  cm  i.d.)  and  were 
extruded  and  frozen  immediately. 

34.  Plant  growth  indices  (stem  biomass,  height,  density  per  m ) 
and  plant-available  nutrient  concentrations  in  sediments  (NH  , PO  , Mn, 
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Ca,  and  Mg)  were  measured  at  all  four  sampling  periods  mentioned  above. 
Two  environmental  factors  were  determined  for  a single  sampling  period 
only — organic  matter  content  of  sediments  (21  July)  and  interstitial 
salinity  (13  September). 

35.  Relative  elevations  of  all  sample  sites  were  determined  with 
a transit,  plane  table,  and  stadia  rod  on  2 July  1975.  Elevations  of 
Taskinas  Marsh  were  rechecked  on  9 September.  Because  no  U.S.G.S.  ele- 
vation bench  marks  were  available  for  reference,  a Mean  Low  Water  datum 
for  both  creeks  was  established  from  tide  height  estimates  published  by 
the  National  Ocean  Survey  (U.S.  Dept,  of  Commerce,  1974).  The  eleva- 
tion of  the  creek  surface  at  slack  low  water  was  measured  relative  to 
the  sample  zones,  then  corrected  for  departure  from  true  mean  low 
water  as  gauged  at  nearby  NOAA  tide  stations.  Slack  low  water  level 
was  surveyed  at  both  creeks  on  14  November  1975,  a day  of  minimal  wind 
activity,  and  rechecked  during  similar  wind  conditions  at  Taskinas 
Creek  on  28  February  1976. 


Laboratory  Procedures 


Plant  samples 


36.  All  live  stems  were  separated  by  species  and  all  dead  stems 
were  also  grouped  together.  The  number  of  live  stems  of  each  species 


was  counted,  and  stem  height  was  measured  to  the  nearest  centimeter 


for  subsamples  of  ten  stems  per  species.  Live  standing  crop  data  for 
each  species  and  total  standing  dead  were  obtained  by  drying  grass 
samples  in  a drying  oven  for  48  hours  at  105°C.  Dried  samples  were 


weighed  to  ± 0.1  gram  on  a Mettler  model  P1200N  top-loading  balance. 
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Sediment  cores 

37.  Analyses  at  the  intensive  sites  were  confined  to  the  upper 
20  cm  of  the  sediment  since  none  of  the  major  plant  species  present 
extends  its  roots  significantly  below  that  level.  (Valiela  and  Teal, 
1973;  Broome  et  al.,  1973).  Soil  subsamples  were  combined  with  a weak 
acid  extractant  solution  (Nelson  et  al.,  1953),  homogenized  in  a 
blender,  and  filtered  to  obtain  samples  for  nutrient  analyses.  Addi- 
tional subsamples  were  ashed  in  a muffle  furnace  to  determine  organic 

content  of  soils  (see  Appendix  A for  details  of  procedures) . 

Nutrient  analysis 

38.  Concentrations  of  iron,  manganese,  calcium,  and  magnesium 

in  filtered  samples  were  measured  in  a Varian  Atomic  Absorption  Spectro- 
photometer (model  AA-5R) . Sample  intake  was  regulated  by  a Technicon 
Sampler  II  automatic  sampling  apparatus,  and  data  were  recorded  as  peak 
heights  on  a Varian  model  A-25  strip  chart  recorder.  Calcium  and  mag- 
nesium were  run  at  60  samples  per  hour  using  the  Technicon  proportion- 


ing pump  to  add  and  mix  one  ml  of  SrCl  (5,000  ppm)  to  each  ml  of  mag- 
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nesium  and  calcium  sample. 

39.  Phosphate  analyses  were  conducted  on  a Technicon  Auto- 
Analyzer;  a combined  reagent  (0.18  M ammonium  molybdate  + .011  M ammon- 
ium metavandate)  was  used  (adapted  from  Steckel  and  Flannery,  1966). 

The  Sampler  was  set  at  40  samples  per  hour,  and  peak  heights  recorded 
on  a Technicon  "Recording  System"  strip  chart.  Duplicate  determina- 
tions were  performed  on  all  the  foregoing  automated  analyses  of  each 
sample. 
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40.  Ammonia  concentrations  were  measured  with  an  Orion  model  95- 
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10  ammonia  electrode  and  an  Orion  model  401  specific  ion  meter.  Two  ml 
of  10  M NaOH  were  added  to  20  ml  of  constantly  stirred  sample  solution 
to  obtain  a pH  of  approximately  11.00;  after  equilibration,  ammonia  con- 
centration was  measured  to  ± 0.1  ppm. 
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PART  V:  FIELD  STUDY  RESULTS  AT  INTENSIVE  SITES 

Habitat  Factors 

Vegetation 

41.  Four  species  of  macrophytic  phanerogams  were  encountered  at 
both  study  sites:  j>.  altemif lora,  S.  patens , D.  spicata,  and  Aster 
tenuifolius  (Figure  4) . j5 . alterniflora  predominated  in  the  low  marsh 
zone  fringing  the  tidal  creeks,  and  j5.  patens  and  D.  spicata  occupied 
adjacent  high  marsh  areas.  All  chree  species  inhabited  the  ecotone  sep- 
arating the  zones,  in  varying  proportions.  A.  tenuifolius,  a minor 
constituent  throughout,  occurred  most  frequently  in  the  low  marsh  and 
ecotone. 

42.  Seasonal  data  for  live  standing  crop,  standing  dead,  stem 
density,  and  stem  height  of  all  species  were  condensed  into  annual  sum- 
mary values  (Figures  5-9).  Because  of  its  relative  scarcity,  A.  tenui- 
folius was  omitted  from  the  summaries.  Maximum  live  standing  crop  of 
the  three  major  species  was  graphed  for  Taskinas  Marsh  (Figure  5a)  and 
Bennett's  Marsh  (Figure  5b).  At  both  sites,  standing  crop  of  j>.  altern- 
iflora generally  declined  from  the  low  marsh  toward  the  high  marsh 
(left  to  right  on  the  graphs) . Maximum  values  occurred  in  the  sample 
zones  closest  to  the  creeks.  j>.  patens  live  standing  crop  exceeded 
D.  spicata  over  most  of  their  common  range  at  both  sites,  although  j>. 
patens  decreased  both  landward  and  seaward  from  peak  values.  Miller 
and  Egler  (1950)  identified  the  J5 . patens  - D.  spicata  association  in 
Connecticut  tidal  marshes;  they  termed  the  species  "ecological 
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STATION  NUMBER 

b.  BENNETT'S  MARSH 

Figure  5.  Maximum  values  for  live  standing  crop  at  Taskinas  Marsh 
(upper)  and  Bennett's  Marsh  (lower). 


equivalents",  while  affirming  that  j3.  patens  is  usually  dominant  where 
they  occur  together.  However,  D.  spicata  surpassed  j>.  patens  along  the 
seaward  margin  of  their  common  range  at  Taskinas;  indeed,  at  Bennett's 
Marsh,  J5.  patens  was  entirely  absent  from  sample  zone  5,  the  seaward  limit 
of  I).  spicata  growth  (Figure  5)  . 

43.  It  is  noteworthy  that  the  largest  standing  crop  values  of  S_. 
patens  exceeded  those  of  j>.  altemif lora  at  both  marshes — a reversal  of 
the  pattern  reported  from  many  other  Virginia  marshes  (Keefe,  1972). 

44.  The  spatial  distribution  of  total  standing  dead  material 
(Figure  6)  is  interesting  for  several  reasons: 

a)  Accumulation  of  standing  dead  is  considerably  greater  in 
the  high  marsh  zone  (S.  patens,  D.  spicata)  than  in  the  low  marsh  (S. 
altemif  lora) . 

b)  The  landward  edge  of  the  ecotone,  adjoining  the  high 
marsh  community,  is  marked  by  an  abrupt  change  in  the  amount  of  stand- 
ing dead. 

c)  The  ratio  of  live  standing  crop  to  dead  standing  crop 
(L/D)  in  both  low  marsh  and  high  marsh  zones  was  considerably  higher 

at  Bennett's  than  at  Taskinas.  Since  L/D  ratios  are  considered  to  vary 
directly  with  the  amount  of  tidal  inundation  that  a marsh  receives 
(Keefe  and  Boynton,  1973;  Mendelssohn,  1973),  Bennett's  Marsh  evidently 
has  experienced  more  inundation  than  Taskinas  Marsh.  Predicted  inunda- 
tion values  (discussed  below)  bolster  this  contention. 

45.  Data  from  July  and  September  sampling  periods  (around  the 
end  of  the  growing  season  in  Virginia  marshes)  were  combined  to  yield 
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average  maximum  values  for  stem  density  and  stem  height.  These 
averages  were  adopted  to  approximate  the  maximum  standing  crop  of  the 
plants  and  to  reduce  the  amount  of  variability  in  the  data. 

46.  The  trend  of  spatial  variation  in  stem  density  was  similar 

in  both  marshes  (Figures  7 and  8) ; stem  density  reached  a maximum  in  the 
middle  of  the  S.  alterniflora  range,  followed  by  a precipitous  decrease 
at  the  ecotone.  Comparison  with  Figure  5 indicated  that  maximum  stand- 
ing crop  was  not  strictly  coincident  with  maximum  stem  density,  although 
their  patterns  of  fluctuation  were  qualitatively  similar.  Stem  densi- 
ties for  the  j>.  patens  - D.  spicata  association  (Figure  8)  also 
resembled  the  respective  standing  crop  trends  (Figure  5).  in  both 
marshes,  D.  spicata  exhibited  its  highest  standing  crop  and  stem  density 
values  at  the  landward  extremes  of  the  transects. 

47.  Stem  heights  (Figure  9),  like  stem  densities,  represented 

averages  of  July  and  September  data.  S^.  alterniflora  steadily  declined 

in  height  rearward;  the  zone  closest  to  the  creek  possessed  both  the 

tallest  stems  and  the  largest  standing  crop  values.  While  j>.  patens 

and  D.  spicata  were  similar  in  stem  height,  j>.  patens  was  slightly 

taller  over  most  of  their  common  range.  Following  the  trend  for 

standing  crop  and  stem  density,  D.  spicata  was  taller  than  j>.  patens  at 

the  seaward  margin  of  the  ecotone  at  Taskinas  Marsh. 

Elevation 

48.  Profiles  of  average  elevation  at  Taskinas  Marsh  and  Bennett's 
Marsh  (Figure  10),  following  the  transects  from  low  marsh  (on  the  left) 
to  high  marsh  (on  the  right),  offer  several  contrasts: 


39 


* STEMS  / M2 


STATION  NUMBER 
Tasklnas  Marsh 


Figure  7.  Maximum  values  for  stem  density  of  Spartina  alternlflora 
at  Taskinas  and  Bennett's  Marsh  (sheet  1 of  2). 
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b.  Bennett9 8 Marsh 
Figure  7 (sheet  2 of  2) 


9 


100B0.0 


S.  PATENS 

B..-SEKAIA 


/\ 


/ 

''o. 


STATION  NUMBER 
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Figure  8.  Maximum  values  for  stem  density  of  Spartina  patens  and 
Distichlis  spicata  at  Taskinas  Marsh  and  Bennett's 
Marsh  (sheet  1 of  2). 
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Figure  8 (sheet  2 of  2), 
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Figure  9 (sheet  2 of  2). 
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Figure  10.  Elevation  of  the  sample  stations  above  mean  low  water. 


a)  Elevation  at  Taskinas  Marsh  increased  sharply  rearward 
from  the  low  marsh  (zones  1-4) , then  leveled  off  in  the  high  marsh 
(zones  5-8).  Bennett's  Marsh  also  rose  toward  the  rear,  but  with  a more 
gradual,  constant  slope  than  Taskinas  Marsh. 

b)  The  slightly  mounded  profile  of  Bennett's  low  marsh,  cen- 
tered on  zone  3,  was  a natural  levee.  Levees  are  common  features  of 
sloping  marsh  foreshores  (Kurz  and  Wagner,  1957).  No  such  structure 
existed  in  the  Taskinas  low  marsh. 

c)  Most  of  the  Taskinas  site  was  considerably  higher  above 

mean  low  water  than  the  Bennett's  Creek  site. 

Number  of  inundations 

49.  A FORTRAN  computer  program  (Program  TIDE)  translated  eleva- 
tion data  into  estimates  of  the  number  of  inundation  events  per  year  at 
each  sample  zone.  An  inundation  event  was  defined  as  any  instance  in 
which  tidal  water  covered  the  sediment  surface.  The  results  can  be  sum- 
marized as  follows  (Figure  11) : 

a)  The  sharpest  deline  in  number  of  inundations  occurred 
between  zones  1 and  2 (the  zones  closest  to  the  tidal  creek)  in  each 
marsh. 

b)  Inundation  events  declined  rapidly  away  from  the  creek 
at  Taskinas  low  marsh;  in  fact,  no  inundation  at  all  is  predicted  land- 
ward of  zone  3.  At  Bennett's  Marsh,  inundation  dropped  more  gradually 
toward  the  rear,  and  only  the  two  highest  zones  (9  and  10)  were  above 
the  reach  of  predicted  tides. 

50.  Thus,  predicted  tidal  inundation  was  considerably  more 
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Figure  11.  Number  of  inundations  per  year  at  sample  stations. 
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extensive  at  Bennett's  Marsh  than  at  Taskinas  Marsh.  As  previously 

stated,  the  fact  that  L/D  ratios  (living: dead  standing  crop)  at 

Bennett's  exceeded  those  at  Taskinas  in  both  low  marsh  and  high  marsh 

corroborates  the  prediction  of  greater  flooding  at  the  Bennett's  Creek 

site  (Mendelssohn,  1973). 

Salinity 

51.  Interstitial  salinity  (Figure  12)  increased  relatively 
sharply  across  the  low  marsh  at  Taskinas  (zones  1-4),  then  peaked  in  the 
ecotone  (zones  4-5)  and  declined  in  the  high  marsh.  Salinity  at  Bennett's 
Marsh  rose  initially  to  a small  peak  at  zone  2,  then  declined  slightly 
before  increasing  again  in  the  high  marsh. 

52.  Although  they  are  not  shown  by  the  averaged  values  in  Figure 
12,  maximum  soil  salinities  exceeded  20  o/oo  at  both  sites,  a concentra- 
tion exceeding  the  highest  values  measured  in  the  creek  water  (about 
15  o/oo).  Soil  salinities  exceeding  the  maximum  salinity  of  tidal  water 
have  also  been  noted  in  other  brackish  estuarine  marshes  (Lindberg  and 

Harriss,  1973). 

Organic  matter 

53.  Samples  from  both  marshes  showed  a nearly  linear  increase  in 
percent  organic  matter  from  the  low  marsh  through  the  high  marsh  (Figure 
13).  Taskinas  Marsh  had  the  greater  range,  possessing  both  the  lowest 
value  (zone  1)  and  the  highest  value  (zone  8).  The  rearward  increase  in 
organic  matter  at  Bennett's  Marsh  was  more  gradual  and  somwehat  less 
linear  than  at  Taskinas  Marsh.  Buttery  et  al.  (1965)  noted  a similar 
rise  in  organic  matter  content  from  low  marsh  to  high  marsh. 

54.  The  organic  content  of  tidal  marsh  sediments  is  largely  a 
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Figure  13.  Percent  organic  matter  of  marsh  sediments. 
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function  of  the  amount  of  plant  material  produced  ^n  situ,  together 
with  the  effects  of  consumption,  decay,  sedimentation,  and  tidal  flush- 
ing. Assuming  that  consumption  and  decay  rates  were  roughly  comparable 
at  Taskinas  Marsh  and  Bennett's  Marsh,  the  differences  in  organic  con- 
tent may  be  ascribed  to  variations  in  productivity,  sedimentation,  and 
tidal  influences.  Shoot  biomass  in  the  high  marsh  community  consid- 
erably exceeded  that  of  the  low  marsh  (Figures  5a  and  5b).  The  fre- 
quency of  tidal  incursions  (Figure  11)  and,  concomitantly,  the  rate  of 
sediment  deposition  (Broome  et  al. , 1973)  were  reduced  toward  the  rear 
marsh.  Both  these  factors  may  help  account  for  the  landward  rise  in 
organic  content.  The  fact  that  Taskinas  high  marsh  produced  more  plant 
tissue  than  Bennett's  high  marsh  and  received  less  tidal  flooding  rela- 
tive to  its  low  marsh  may  explain  the  steeper  front-to-rear  rise  in 
organic  matter  content  at  Taskinas. 


Plant-Available  Nutrients 

rv 

55 . Summary  graphs  presented  and  discussed  in  the  following 

. 


paragraphs  represent  averages  of  three  transects,  all  four  sampling 


times,  and  both  depth  ranges  (0-10  cm  and  10-20  cm);  that  is,  each  data 


point  on  the  figures  encompasses  24  individual  values  obtained  by  sedi- 


ment analysis.  The  means  and  standard  errors  for  the  habitat  variables 

are  given  in  Table  1. 

Ammonia-N 

56.  Summary  values  for  estimated  plant-available  ammonia-N 
(Figure  14)  indicated  that,  at  both  marshes,  slightly  higher  values 
tended  to  characterize  the  high  marsh  sediments.  However,  the  range  of 


Table  1 

Means  (X)  and  Standard  Errors  (S.E.M.)  of  Habitat  Parameters  from  the  Study  Sites 


Note:  Nutrient  data  include  4 sampling  periods  (4-5  April,  29-30  May,  21  July,  and  13  September)  and  2 depth  ranges  in  sediment  (0-10  cm  and  10-20  cm). 
Sample  tones  are  numbered  consecutively  along  transects  from  low  marsh  to  high  marsh;  transect  length  » 28  m (Tasklnas  Marsh)  and  36  m (Bennett’s 


STATION  NUMBER 

Figure  14.  Concentration  (ppm)  of  plant-available  ammonia-N  in 
marsh  sediments  at  sample  stations. 


variation  was  small  (5. 5-9.0  ppm)  and  clear  patterns  of  spatial  change 
were  not  evident. 


57.  Jeffrey  (in  Piggott,  1969)  also  found  that  concentrations  of 
ammonia  exhibited  little  pattern  across  a marsh.  However,  since  inor- 
ganic nitrogen  can  change  quite  rapidly  in  form  and  amount  in  soils  (by 
microbial  uptake  and  release,  plant  absorption,  volatilization,  immobil- 
ization in  humus,  absorption  within  crystal  lattices,  and  mineraliza- 
tion of  organic  matter)  concentrations  of  ammonia-N  at  a few  instants 

of  time  may  not  accurately  represent  the  soil's  capacity  to  supply  inor- 
ganic nitrogen  to  plants  over  an  entire  growing  season  (Tusneem  and 

Patrick,  1971). 

Phosphate 

58.  Both  marshes  displayed  a relatively  rapid  increase  in  avail- 
able phosphate  concentration  from  zone  1 to  zone  2 (Figure  15).  Piggott 
(1969)  also  reported  a landward  increase  in  phosphorus  content  of  marsh 
soils.  At  Taskinas,  the  values  continued  to  increase  rearward  but  at  a 
slower  rate,  leveling  off  in  the  high  marsh.  The  Bennett's  Creek  tran- 
sect showed  two  small  peaks,  one  each  in  the  high  marsh  and  low  marsh, 
with  a slight  trough  toward  the  rear  of  the  ecotone(zone  6).  Concentra- 
tions at  Taskinas  were  generally  somewhat  lower  than  at  Bennett's  Creek. 
Iron 

59.  Concentrations  of  plant-available  iron  declined  gradually 
from  the  low  marsh  toward  the  rear  zones  at  both  study  sites  (Figure 
16).  At  Taskinas,  the  concentrations  were  somewhat  higher  and  the  rear- 
ward decline  was  less  linear  than  at  Bennett's.  However,  zones  9 and 

10  mark  a major  departure:  a sudden  steep  increase  in  concentration, 
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attaining  the  highest  values  recorded  for  Bennett's  Marsh. 

60.  The  results  varied  with  regard  to  Adams'  (1963)  contention 
that  tall  E^.  alternif lora  tended  to  be  associated  with  higher  iron  con- 
centrations than  the  S.  patens  - I).  spicata  association.  The  data  from 
Taskinas  Marsh  agreed  with  Adams'  findings,  but  Bennett's  Marsh  data 
indicated  exactly  the  opposite  pattern:  higher  iron  concentrations  in 
the  high  marsh.  Direct  comparisons  with  Adams  may  be  tenuous,  however; 
he  did  not  specify  the  extraction  procedures  he  used,  and  the  iron  con- 
centrations recorded  from  his  North  Carolina  marshes  were  more  than  an 

order  of  magnitude  lower  than  those  from  the  present  study. 

Manganese 

61.  Both  sites  displayed  a rapid  decrease  in  available  manganese 

concentration  from  the  low  marsh  toward  the  high  marsh,  with  nearly  all 

the  reduction  confined  to  the  first  few  zones  (Figure  17).  At  Taskinas 

Marsh,  however,  the  highest  value  (zone  1)  was  more  than  six  times 

greater  than  the  maximum  at  Bennett's  Marsh,  and  the  initial  decline 

was  much  more  precipitous.  Concentrations  from  the  ecotone  landward 

were  low  and  variability  was  minimal;  averages  at  Bennett's  Marsh 

fluctuated  within  only  0 to  4 ppm  between  zone  4 and  zone  9. 

Calcium 

62.  Calcium  concentrations  at  the  two  marshes  varied  in  a very 
similar  manner  (Figure  18) . Both  exhibited  a peak  value  in  the  low 
marsh  adjacent  to  the  ecotone,  a decline  rearward  through  the  ecotone, 
and  some  upward  recovery  in  the  high  marsh  (the  latter  trend  was  more 
evident  at  Bennett's  Creek  Marsh). 
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STATION  NUMBER 

Figure  17 . Concentration  (ppm)  of  plant-available  manganese  in 
marsh  sediments  at  sample  stations. 


i?3t*«rt«ii 


STATION  NUHBER 

Figure  18.  Concentration  (ppm)  of  plant-available  calcium  in 
marsh  sediments  at  sample  stations. 
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Magnesium 

63.  The  graphs  of  available  magnesium  (Figure  19)  showed  a strik- 
ing resemblance  to  the  pattern  for  calcium  (Figure  18).  As  with  cal- 
cium, both  marshes  varied  in  nearly  identical  fashion:  high  values 
occurred  at  the  seaward  edge  of  the  ecotone,  and  Bennett's  Marsh 
attained  its  highest  concentrations  in  zones  9 and  10. 

Correlation  Analyses 

64.  Statistical  correlation  is  a useful  first  step  toward  identi- 
fying the  habitat  factors  most  closely  related  to  plant  growth  and  dis- 
tribution. Bivariate  correlation  is  a statistical  technique  for  mea- 
suring the  intensity  of  relationship  between  co-occurring  variables. 

The  correlation  statistic  employed  here  is  the  Pearson  product-moment 
correlation  coefficient,  "r". 

65.  The  test  used  assumes  that  the  data  are  randomly  sampled, 
normally  distributed,  and  there  is  a linear  relation  between  the  vari- 
ables. The  data  set  was  sufficient  to  generally  satisfy  the  latter  two 
correlations.  The  data  were  sampled  by  a form  of  stratified  random 
method,  whereby  samples  were  collected  randomly  from  within  specified 
sample  areas.  Values  of  "r"  and  their  associated  statistical  signifi- 
cance were  calculated  with  the  "Scattergram"  subprogram  of  Statistical 
Package  for  the  Social  Sciences  (SPSS)  (Nie  et  al. , 1975). 

66.  Tables  2 and  3 display  correlation  coefficients  relating 
plant  growth  variables  to  habitat  variables  at  Taskinas  Marsh  and 
Bennett's  Marsh,  respectively.  The  three  growth  variables  measured— 
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Figure  19.  Concentration  (ppm)  of  plant-available  magnesium  in 
marsh  sediments  at  sample  stations. 
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standing  crop,  stem  height,  and  stem  density — are  listed  for  each  of 
the  three  predominant  plant  species:  S.  alternif lora,  j>.  patens,  and 
D.  spicata. 

67.  Several  of  the  plant-habitat  correlations  deserved  special 
emphasis.  It  was  evident  from  comparing  spatial  patterns  of  standing 
crop  (Figure  5)  and  stem  height  (Figure  8 ) with  elevation  and  inunda- 
tion (Figures  10  and  11)  that  growth  of  alterniflora  exhibited  a 
highly  significant  positive  correlation  with  tide  factors.  This  rela- 
tionship had  been  observed  repeatedly  (Johnson  and  York,  1915;  Adams, 
1963;  Kurz  and  Wagner,  1957);  but  since  the  sampling  sites  in  the  pre- 
sent study  did  not  extend  to  the  seaward  extreme  of  j5.  alterniflora 
distribution,  the  contention  that  tidal  inundation  restricted  the  spe- 
cies along  the  lower  margin  of  its  range  within  the  marsh  (Johnson  and 
York,  1915;  Teal  and  Kanwisher,  1961)  could  not  be  tested.  For  S^. 
patens  and  D.  spicata,  the  measured  responses  to  tides  differed  conspi- 
cuously between  the  two  marshes.  At  Taskinas,  where  tidal  water  very 
seldom  reached  the  high  marsh  (Figure  10  ) growth  of  neither  species 
was  closely  related  to  inundation  or  elevation;  Bennett's  high  marsh, 
which  received  some  regular  flooding  during  spring  tides,  showed  close 
correlations  of  tide  with  j>.  patens  (positive)  and  D.  spicata  (negative) . 

68.  As  indicated  by  previous  investigators  (Gosselink,  1970; 
Palmisano,  1970;  Phleger,  1971),  growth  of  certain  species  tends  to  be 
inhibited  by  increasing  soil  salinity.  Negative  correlations  were  sig- 
nificant for  j5.  alterniflora  at  both  marshes  and  less  pronounced  for 

S.  patens.  D.  spicata,  the  most  salt-tolerant  of  the  three  species 
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(Adams,  1963),  showed  positive  correlations  with  salinity  at  Bennett's 
Marsh. 

69.  The  absence  of  significant  correlations  between  ammonia-N 
concentration  and  growth  of  any  species  appears  to  support  Tusneem  and 
Patrick's  (1971)  contention  that  sampling  only  at  isolated  instances 
(as  was  done  in  this  study)  may  misrepresent  the  plant  response  to  soil 
nitrogen  over  a growing  season.  The  fact  that  nitrogen  availability 
clearly  limits  growth  of  marsh  grasses  (Valiela  and  Teal,  1973;  D.E. 
Marshall,  1970;  Broome  et  al.,  1973),  coupled  with  the  low  concentra- 
tions found  in  the  soil  (Figure  14).  suggested  that  rapid  depletion  by 
plant  uptake  may  engender  considerable  fluctuation  of  ammonia-N  levels. 
Analysis  of  available  nitrogen  in  both  soil  and  plant  tissue  at  fre- 
quent intervals  over  the  growing  season  might  provide  a better  measure 
of  plant  response  to  soil  nitrogen.  Working  with  aquatic  species, 
Gerloff  and  Krombholz  (1966)  noted  a critical  level  of  nitrogen  concen- 
tration in  plant  tissues  of  1.3  percent  below  which  growth  was  inhib- 
ited; similar  data  for  salt  marsh  species  would  help  elucidate  the  rela- 
tionship between  plant  growth  and  nitrogen  availability  in  soils. 

70.  Phosphate  concentrations  were  not  closely  correlated  with 
plant  growth,  either.  Since  fine-textured,  highly  reduced  sediments, 
such  as  typify  Taskinas  Marsh  and  Bennett's  Marsh,  usually  possess  abun- 
dant available  phosphorus  (Broome  et  al. , 1973),  the  absence  of  a close 
relationship  with  plant  growth  was  not  surprising. 

71.  j>.  altemlf lora  standing  crop  correlated  significantly  with 
iron  availability  at  both  marshes,  supporting  the  theory  that  this 


species  thrives  where  iron  is  plentiful  (Adams,  1963;  Mooring  et  al., 
1971;  Nixon  and  Oviatt,  1973).  JS . patens  and  D.  spicata  showed  opposite 
patterns  of  correlation  at  Bennett’s  Marsh:  S.  patens  was  negatively 
correlated  to  iron,  while  D.  spicata  had  a highly  significant  positive 
correlation. 

72.  A contrast  is  discernible  between  Table  2 and  Table  3 in  the 
arrangement  of  significant  correlations.  Significant  plant-habitat  cor- 
relations were  almost  entirely  restricted  to  j>.  alterniflora  at 
Taskinas  Marsh,  while  all  three  species  exhibited  significant  correla- 
tions at  Bennett's  Marsh.  This  difference  suggests  that  more  of  the 
habitat  factors  measured  were  closely  related  to  growth  of  high  marsh 
species  (S.  patens  and  I).  spicata)  at  Bennett's  Marsh  than  at  Taskinas 
Marsh. 

73.  Figure  10  and  Figure  11  indicate  that  while  Taskinas  Marsh 
received  no  regular  tidal  influence,  Bennett's  high  marsh  was  inundated 
to  some  extent.  This  pattern  of  correlations  appeared  to  support  the 
hypothesis  (Hinde,  1954;  Gray  and  Bunce,  1972)  that  tidal  changes  in 
water  level  may  exert  a controlling  influence  upon  plant  growth  by  modi- 
fication of  other  habitat  conditions. 

74.  Tables  4 and  5 are  matrices  of  correlation  coefficients  that 
indicate  the  intensity  of  relationships  among  habitat  factors.  Values 
of  "r"  and  their  significance  were  calculated  for  the  distinct  communi- 
ties in  each  marsh.  Low  marsh  correlations  included  only  the  range  of 
j3.  alterniflora,  while  high  marsh  correlations  were  restricted  to  the 
J3.  patens  - D.  spicata  community. 


* ” significant  at  5Z  level;  **  - significant  at  1Z  level. 

Note:  Values  above  the  main  diagonal  are  for  the  high  marsh;  values  below  the  main  diagonal  are  for  the  low  marsh. 


Correlation  Coefficients  between  Habitat  Variables  at  Bennett' s Marsh 


* “ significant  at  5%  level;  **  - significant  at  12  level. 

Note:  Values  above  the  main  diagonal  are  for  the  high  marsh;  values  below  the  main  diagonal  are  for  the  low  marsh. 
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75.  As  a source  of  inorganic  nutrients,  a provider  of  abundant 
sites  for  cation  sorption,  and  a fuel  for  microbially  mediated  reac- 
tions, organic  matter  is  a vital  constituent  of  the  soil  environment 
(Buckman  and  Brady,  1969).  Because  reduction  of  manganese  and  iron  to 
plant-available  forms  is  facilitated  by  anaerobic  bacteria  (Turner  and 
Patrick,  1968),  strong  positive  correlations  between  organic  content  and 
these  cations  might  be  exnected.  In  fact,  most  of  the  correlations 
were  highly  negative  (Tables  4 and  5).  Nixon  and  Oviatt  (1973)  noted  a 
close  positive  correspondence  between  iron  concentration  and  percent 
organic  matter,  but  the  organic  content  of  their  Connecticut  marshes 
averaged  less  than  half  that  found  in  the  present  research.  Conceivably, 
available  cations  in  the  highly  organic  soils  of  Taskinas  Marsh  and 
Bennett's  Marsh  were  rapidly  incorporated  into  microbial  biomass,  and 
thus  became  unavailable  to  plants.  Organic  matter  was  not  the  sole 

influence  on  iron  and  manganese  availability;  however,  long  periods  of 

2+  2*f‘ 

waterlogging  tended  to  liberate  large  amounts  of  available  Mn  and  Fe 
(Leeper,  1947);  indeed,  a strong  positive  correlation  existed  between 
tidal  inundation  and  levels  of  these  cations,  particularly  in  the  low 
marsh  zone  (Tables  4 and  5). 

76.  Many  species  of  cations  are  influenced  by  the  behavior  of  iron 
in  marsh  soils.  When  surface  sediments  become  aerated  during  tidal 
emergence,  a number  of  trace  elements,  including  manganese, 
co-precipitate  with  iron  in  the  presence  of  organic  matter 
(Bloomfield,  1963).  Like  manganese,  iron  is  readily  reduced  by  microbes. 
The  consistently  high  positive  correlations  between  iron  and  manganese 
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(Tables  4 and  5)  further  indicate  that  the  behavior  of  iron  and  mangan- 
ese in  marsh  soils  was  closely  allied.  Iron  is  also  directly  involved 
in  the  absorption  of  phosphorus  by  marsh  soils  (Jitts,  1959).  Hydrated 
iron  oxides  provide  binding  sites  for  phosphorus  ions;  ferrous  phosphate 
is  believed  to  be  an  important  mobile  compound  in  reduced  soils 
(Patrick  and  Khalid,  1974).  Iron  and  phosphate  concentrations  showed 
close  positive  correlation  at  Taskinas  low  marsh  and  Bennett's  high 
marsh,  but  elsewhere  the  relationship  was  either  weak  or  strongly  nega- 
tive. Calcium  phosphates  are  also  said  to  be  an  important  phosphorus 
source  in  marsh  soils  (Broome  et  al. , 1973) , but  calcium  and  phosphate 
correlations  were  inconsistent. 

77.  Calcium  and  magnesium  concentrations  were  highly  correlated 
at  both  marshes,  reflecting  the  very  close  resemblance  of  their  spatial 
profiles  (Figures  18  and  19).  Laboratory  studies  indicated  that  both 
cations  increase  in  concentration  during  prolonged  submergence,  and 
that  their  availability  is  similarly  influenced  by  the  decay  of  organic 
matter  (Robinson,  1930). 

78.  The  foregoing  discussion  of  edaphic  interactions  raises  more 
questions  than  it  answers.  However,  it  demonstrates  that  the  interac- 
tions among  sediment  variables  possess  a high  degree  of  complexity  and 
variability — a variability  that  contributes  heavily  to  the  heterogeneity 
of  marsh  soils  and  of  their  associated  vegetation  (Gray  and  Bunce,  1972). 

79.  Further  examination  of  the  habitat  intercorrelations  (Tables 
4 and  5)  reveals  that  the  tidal  factors  (elevation,  inundations  per 
year,  and  hours  of  inundation)  displayed  numerous  significant 
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correlations  throughout  Bennett's  Marsh  and  at  Taskinas  low  marsh,  but 
lacked  significant  correlations  at  Taskinas  high  marsh — the  same  pattern 
that  typified  the  plant-habitat  correlations  (Tables  2 and  3).  Of  all 
the  habitat  factors  considered,  only  the  tide  variables  possessed  sim- 
ilar patterns  of  correlation  with  both  plant  factors  and  other  habitat 
factors. 

80.  In  regularly  inundated  areas,  both  tide  vs.  habitat  and  plant 
growth  vs.  habitat  correlations  were  strong;  above  the  reach  of  the  tides 
plant  growth  is  no  longer  closely  related  tc  the  habitat  conditions  mea- 
sured. This  contrast  lends  circumstantial  support  to  the  theory  that 
tides  act  to  influence  plant  growth  by  controlling  other  aspects  of  the 
marsh  habitat.  However,  it  does  not  diminish  the  likelihood  that  inunda- 
tion also  influences  plant  growth  directly  (Hinde,  1954)  by  imposing 
restrictions  upon  root  respiration,  nutrient  uptake,  seed  germination,  or 
other  biotic  processes.  It  seems  likely  that  direct  and  indirect  tidal 
influences  are  both  important  controls  on  plant  growth. 

Competition  Analyses 

81.  The  present  data  indicate  that  plant  growth  is  not  closely 
related  to  edaphic  conditions  in  the  absence  of  tidal  inundation.  Sev- 
eral investigators  (Johnson  and  York,  1915;  Chapman,  1938;  Reed,  1947) 
have  suggested  that  competition  may  affect  zonation  of  marsh  plants  above 
the  tide  line.  Johnson  and  York  (1915)  theorized  that  the  landward  limit 
of  occurrence  of  J3.  alternif lora  is  determined  primarily  by  competition 
with  species,  such  as  j>.  patens  and  D.  spicata,  whose  inability  to 
tolerate  prolonged  submergence  precludes  their  seaward  expansion.  Reed 
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(1947)  amplified  the  competition  hypothesis  by  suggesting  that  the 
accretion  of  sediments  on  top  of  j>.  alterniflora  at  the  ecotone  might 
confer  an  advantage  upon  high  marsh  species  in  competing  for  growing 
space.  Additional  mechanisms  of  competition  may  be  physical,  such  as 
inhibition  of  photosynthesis  by  shading  of  leaf  surfaces  (Ranwell,  1964; 
Buttery,  1965),  or  chemical,  such  as  allelopathic  interactions.  Examina- 
tion of  Figures  5 and  9 reveals  that  j3.  patens  usually  enjoys  an 
advantage  in  the  competition  for  sunlight. 

82.  Competitive  success  or  failure  is  closely  tied  to  habitat 
phenomena,  yet  even  less  is  known  about  plant  interactions  than  about 


the  nature  of  the  marsh  habitat.  In  fact,  a thorough  quantitative  study 
of  plant  competition  in  marshes  has  never  been  accomplished. 

Gradient  Analysis 

83.  Figures  20  and  21  are  examples  of  community  gradient  graphs 
of  the  study  sites.  On  each  figure,  values  of  a single  habitat  factor 
were  matched  with  corresponding  values  of  maximum  live  standing  crop 
data  for  each  species;  the  standing  crop  data  were  expressed  as  percent- 
ages of  the  largest  value.  Lines  were  drawn  to  indicate  patterns  of 
change  in  standing  crop  (y-axis)  over  the  range  of  the  particular  habi- 
tat factor  (x-axis).  Each  line  represents  a two-point  running  mean  of 
the  data  points. 

84.  The  extent  to  which  a given  habitat  factor  influences  the 
segregation  of  the  two  communities  can  be  estimated  by  the  degree  of 
isolation  of  the  growth  maximum  of  j>.  alterniflora  (low  marsh)  from  the 

growth  maxima  of  _S.  patens  and  D.  spicata  along  the  habitat  axis.  1 
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Figure  20.  Community  gradient  plots  of  live  standing  crop  vs.  eleva- 
tion at  v'a)  Taskinas  Marsh  and  (b)  Bennett's  Marsh 
(sheet  1 of  2) . 
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Figure  21.  Community  gradient  plots  of  live  standing  crop  vs.  salin- 
ity at  (a)  Taskinas  Marsh  and  (b)  Bennett's  Marsh 
(sheet  1 of  2) . 
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Because  the  scales  of  the  habitat  factors  have  not  been  normalized. 


this  evaluation  is  somewhat  arbitrary;  nevertheless,  the  patterns 
represented  can  be  compared  qualitatively.  For  example.  Figure  20 
shows  clearcut  separation  of  the  S^.  alternif lora  response  peak  from  the 
peaks  of  j3.  patens  and  D.  spicata.  In  contrast,  Figure  21  displays 
multiple  response  peaks  for  each  species,  with  poorly  defined  separa- 
tion of  the  j5.  alternif lora  peak  from  the  peaks  of  the  high  marsh 
species.  Comparing  these  results  yields  the  conclusion  that,  in  both 
marshes,  elevation  is  a relatively  more  active  agent  than  ammonia  con- 
centration (not  shown)  in  separating  high  marsh  from  low  marsh  com- 
munities. While  only  Figures  20  (elevation)  and  21  (salinity)  are  pre- 
sented here  for  purposes  of  illustration,  community  gradient  graphs 
were  prepared  for  each  habitat  variable  (tidal  factors  were  represented 
by  the  elevation  variable  only) . 

85.  The  following  variables  manifested  the  clearest  separation  of 


S.  alternif lora  peaks  from  S.  patens  and  D.  spicata  peaks: 


TASKINAS  MARSH 


BENNETT'S  MARSH 


elevation 


elevation 


organic  matter 


organic  matter 


phosphate 


calcium 


manganese 


magnesium 


manganese 


The  community  gradient  technique  thus  singled  out  the  above  factors  as 
the  most  influential  in  dividing  the  vegetational  zones. 


86.  A second  phase  of  gradient  analysis,  the  complex  gradient 
approach  (Whittaker,  1967),  arrayed  pairs  of  these  selected  habitat 
factors  on  opposing  axes  and  plots  corresponding  standing  crop  data  in 
habitat  spaces.  Standing  crop  data  were  plotted  on  a two-level  rating 
scale: 

optimal  - 50  percent  or  more  of  maximum  live  standing  crop 
sub-optimal  - less  than  50  percent  of  maximum  live  standing 
crop 

This  two-part  scale  superseded  the  "percent"  scale  used  on  community 
gradient  graphics. 

87.  The  object  of  the  complex  gradient  analysis  was  to  pinpoint 
pairings  of  habitat  factors  that  most  clearly  separate  the  optimal 
ranges  of  J5.  patens  and  I).  spicata  from  that  of  j>.  alterniflora. 

88.  Illustrative  examples  of  complex  gradient  plots  are  presented 
for  Taskinas  Marsh  (Figures  22  and  23)  and  Bennett' s Marsh  (Figures  22 
and  24).  In  each  plot,  the  optimal  portion  of  a species'  range  is 
shaded  by  parallel  lines,  while  the  sub-optimal  range  is  merely  out- 
lined. 

89.  At  Taskinas  Marsh,  the  combination  of  elevation  and  organic 
matter  (Figure  22)  appeared  to  produce  the  clearest  separation  of  high 
marsh  and  low  marsh  optima  dong  both  axes.  Other  combinations  of 
selected  factors,  such  as  manganese  and  phosphate  (Figure  23),  show  a 
greater  degree  of  overlapping  of  optima.  Similarly,  only  the  graph  of 
elevation  and  organic  matter  at  Bennett's  Marsh  (Figure  22  ) clearly 
separated  the  j>.  alterniflora  optimum  from  the  optima  of  the  high 
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marsh  species  along  both  axes.  The  consistency  of  the  results  between 
the  plot  of  manganese  and  calcium  at  Bennett's  Marsh  (Figure  24)  exem- 
plified factor  combinations  with  overlapping  low  marsh  and  high  marsh 
optima  and  led  to  the  conclusion  that , among  the  habitat  factors  con- 
sidered in  this  study,  elevation  and  organic  matter  were  most  closely 
associated  with  community  zonation. 

90.  Both  organic  matter  content  and  tidal  variables  were  closely 
associated  with  plant  growth  and  distribution  at  the  study  sites;  how- 
ever, their  influences  were  dissimilar  in  nature.  Highly  organic  soils 
often  suppport  very  productive  flora  because  detrital  decay  releases 
inorganic  nutrients  and  abundant  colloidal  particles  bind  exchangeable 
cations  (Albrecht,  1941;  Buckman  and  Brady,  1969).  In  addition,  decay 
of  organic  matter  tends  to  lower  soil  pH,  altering  the  availability  of 
many  plant  nutrients.  Thus,  organic  matter  seems  most  likely  to  affect 
plant  growth  through  its  influence  on  the  availability  of  inorganic 
plant  nutrients  (Bloomfield,  1963). 

91.  Where  plant  growth  is  correlated  to  some  extent  with  avail- 
able nutrient  levels,  one  might  expect  to  find  a similar  degree  of  cor- 
relation between  organic  matter  content  and  these  nutrient  concentra- 
tions. Tables  2 and  3,  together  with  Tables  4 and  5,  indicate  that 
this  correspondence  was  found  in  the  low  marsh  zones,  but  not  in  the 
high  marsh  zones.  Considering  all  three  major  plant  species  together, 
then,  the  influence  exerted  by  organic  matter  upon  plant  growth  appears 
to  be  somewhat  inconsistent.  As  noted  previously,  tide-habitat  corre- 
lation patterns  are  paralleled  by  tide-plant  growth  correlation 
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patterns.  Comparing  tide  and  organic  matter  correlation  patterns  in 
this  manner  suggests  that  tidal  influences  may  be  the  more  pervasive. 
Furthermore,  since  sedimentation  rate,  a determinant  of  the  organic 
matter  content  of  marsh  soils,  is  itself  largely  a function  of  the 
degree  of  tidal  inundation  (Broome  et  al.,  1973),  it  seems  apparent 
that  tides  exert  some  control  over  organic  matter  content.  The  sum  of 
the  evidence  in  this  project,  therefore,  appears  to  point  toward 
tidal  inundation  as  the  major  controlling  influence  upon  plant  growth 
and  zonation. 

92.  Prevalence  of  the  close  relationship  between  tide  factors 
and  vegetation  tends  to  obscure  major  differences  between  the  marshes. 
Certain  habitat  factors  were  significantly  correlated  with  plant 
growth  at  only  one  marsh  (Tables  2 and  3);  interrelationships  among 
habitat  factors  (Tables  4 and  5)  also  varied  between  Taskinas  Marsh 
and  Bennett's  Marsh. 

93.  Even  the  zonations  of  species  with  regard  to  tidal  inunda- 

* tion  were  somewhat  dissimilar.  Table  6 contrasts  the  two  marshes  in 

terms  of  species  tolerance  to  tidal  inundation.  The  two  inundation 
parameters  were  calculated  only  for  the  extremes  of  species  distribu- 
tion that  fell  within  the  sample  zone;  that  is,  the  landward  limit  of 
j>.  alterniflora  and  the  seaward  limits  of  j>.  patens  and  I).  spicata. 

O'i . Inundation  tolerances  were  uniformly  greater  at  Bennett's 
Marsh  than  at  Taskinas  Marsh.  j>.  alterniflora  tolerated  a minimum  of 
six  inundations  per  year  at  Bennett's,  but  survived  with  no  inundation 
at  Taskinas.  Both  patens  and  D.  spicata  survived  considerably 
more  inundation  at  the  lower  limits  of  their  range  at  Bennett's  than 
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Table  6 

Limits  of  Tolerance  to  Tidal  Inundation  for  Three  Species  of 
Marsh  Grasses  at  Taskinas  Marsh  and  Bennett ' s Marsh 


Inundation  Events  per  Yr.  Hrs.  of  Inundation  per  Yr. 


Species 

Taskinas 

Bennett's 

Taskinas 

Bennett's 

Max. 

Min. 

Max . Min . 

Max.  Min. 

Max.  Min. 

S.  alterniflora 

— 

0 

6 

0 

3.5 

S.  patens 

8 

0 

17 

6.5 

22 

D.  splcata 

8 

— 

41 

6.5 

65 

Note:  Max.  = maximum  number  of  events  recorded 
Min.  = minimum  number  of  events  recorded 
* not  measured  within  sampling  sites 


86 


at  Taskinas.  While  D.  spicata  tolerated  several  times  more  inundation 
than  S^.  patens  at  Bennett's,  the  two  species  showed  the  same  tolerance 
limits  at  Taskinas. 

95.  It  should  be  noted,  however,  that  annual  estimates  of  tidal 
inundation  are  subject  to  inaccuracy  from  several  environmental  sources, 
including  alteration  of  the  tidal  prism  by  wind,  hindrance  of  tidal 
waters  by  physical  resistance  of  the  vegetation,  and  variability  of  sed- 
iment drainage  rates  (Kurz  and  Wagner,  1957).  Inundation  of  Taskinas 
high  marsh,  though  not  predicted,  has  been  observed  when  high  tides 

coincide  with  strong  northwest  winds  (J.  Pickral,  1976,  personal  com- 
* 

munication  ) . 

96.  Despite  estimation  errors,  tidal  tolerances  of  each  plant 
species  appeared  to  differ  markedly  between  the  two  marshes.  Assum- 
ing that  a given  species  has  the  genetic  potential  for  response  to 
tidal  inundation  in  both  marshes  (Mooring  et  al.,  1971),  differences  in 
response  might  be  explained  by  the  divergent  past  histories  of  the 
marshes.  As  discussed  previously,  Taskinas  Marsh  is  presently  under- 
going a "drowning"  phase,  while  Bennett's  Marsh  appears  to  be  building 
upward  by  sediment  accretion.  A gradually  sloping  foreshore,  such  as 
at  Bennett's  Marsh,  is  frequently  associated  with  a young,  expanding 
marsh  (Redfield,  1972),  while  the  steeper  foreshore  at  Taskinas  is  more 
typical  of  a mature  marsh  (Stevenson  and  Emery,  1958).  The  difference 
profoundly  affects  the  extent  of  tidal  inundation  (Figure  11).  The 

*J.  Pickral,  1976,  Dept  of  Environmental  Sciences,  University  of 
Virginia,  Charlottesville,  VA  22903. 
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absence  of  predicted  inundation  in  Taskinas  high  marsh  may  result  from 
a relatively  recent  high  stand  of  the  marsh,  as  evidenced  by  the  pre- 
sence of  scattered  dead  _I.  frutescens  among  the  high  marsh  grass.  In 
other  words,  the  contemporary  tidal  regime  of  the  high  marsh  was  condi- 
tioned by  this  past  episode  as  well  as  by  the  current  drowning  phase. 

97.  In  mature  marshes,  particularly,  vegetational  succession  in 
response  to  changing  relative  sea  level  experiences  time  lags.  Reed 
(1947)  suggested  that,  in  an  aggrading  marsh,  j>.  alterr.iflora  stems 
buried  by  accumulating  sediments  may  remain  alive  for  some  time,  since 
their  roots  continue  to  exist  under  relatively  unchanged  conditions. 
Eventually,  the  high  marsh  species  encroaching  on  the  new  sediments 
smother  the  j3.  alterniflora  beneath  them,  and  the  transition  is  com- 
pleted. Similarly,  the  retrogressive  succession  accompanying  rising 
relative  sea  level  (Kurz  and  Wagner,  1957)  often  requires  that  an  invad- 
ing species  (e.g.,  S.  alterniflora  at  Taskinas  Marsh)  establish  itself 
among  a dense,  fibrous  mat  of  high  marsh  vegetation.  Ranwell  (1954) 
observed  that  competition  from  Spartina  townsendii  appeared  to  retard 
the  advance  of  invading  species,  even  when  habitat  conditions  would 
have  permitted  the  invaders  to  expand  more  rapidly.  Pre-existing  vege- 
tation thus  may  possess  a competitive  advantage  by  virtue  of  its  physi- 
cal occupation  of  the  growing  space,  enabling  it  to  persist  for  some 
time  despite  an  increasingly  unfavorable  habitat.  The  equilibration 
time  of  marsh  vegetation  during  periods  of  change  in  relative  sea  level 
may  significantly  modify  the  relationship  between  vegetational  distri- 
bution and  tidal  inundation. 
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PART  VI:  MONITORING  RESULTS 


98.  As  the  monitoring  equipment  was  built,  it  was  tested  success- 
fully in  the  laboratory;  in  the  harsher  environment  of  the  salt  marsh, 
however,  unforeseen  problems  arose.  Throughout  the  field  monitoring 
phase,  the  equipment  was  plagued  by  component  failure  of  varying 
severity.  The  initial  difficulties  were  most  serious;  failure  of  one 
or  more  components  would  induce  problems  in  others,  and  the  system's 
performance  would  totally  deteriorate.  Chief  among  these  early  prob- 
lems were  loss  of  amplifiers  due  to  electronic  failure,  and  disappear- 
ance of  a transmitted  signal  or  shorting  out  due  to  water  leakage 
around  the  submersible  connectors. 

99.  When  the  amplifiers  malfunction,  the  unit  frequently  sets  up 
an  overvoltage.  The  recorder  normally  receives  an  input  of  0.0  to  1.0 
v-d.c.;  any  higher  voltage  produces  erroneous  records  on  the  channel  in 
which  the  overvoltage  occurs  and  in  adjacent  channels  as  well.  A suffi- 
ciently high  overvoltage  damages  the  recorder  itself.  Ultimately,  volt- 
age regulators  were  installed  on  the  outputs  of  all  amplifiers  to  pre- 
vent overvoltages,  but  many  records  were  lost  before  this  problem  was 
alleviated. 

100.  The  other  major  problem,  which  was  not  fully  resolved, 
involved  the  shorting  out  of  signals  between  the  amplifiers  and  the 
junction  boxes.  The  major  cause  of  this  problem  proved  to  be  leakage 
in  the  submersible  cable  connectors.  Since  it  was  frequently  necessary 
to  carry  the  equipment  several  hundred  meters  to  the  study  site,  the 


units  were  designed  in  a modular  fashion  to  facilitate  handling  assembly 


on  the  marsh.  Despite  careful  handling,  bits  of  sediment  and  droplets 
of  salt  water  frequently  became  attached  to  exposed  connectors  and  were 
not  completely  removed  by  a final  cleaning  prior  to  assembly  of  the  con- 
nectors. Any  trace  of  water  or  salt  on  the  connectors  allowed  the  vol- 
tage to  leak,  either  to  ground  or  to  another  connector.  This  leakage 
periodically  caused  erroneous  behavior  which  was  often  perplexing  and 
difficult  to  trace.  Despite  the  host  of  electronic  problems  in  the 
system  and  the  consequent  lack  of  complete  records,  the  results  from 
the  monitoring  system  were  quite  encouraging.  The  main  value  of  this 
part  of  the  project  was  the  demonstration  of  the  feasibility  of  record- 
ing small  time  interval  variations  in  sedimentary  parameters  and  relat- 
ing these  variations  to  plant  distribution. 

101.  The  most  complete  set  of  data  collected  with  the  system  was 
collected  between  25  and  28  November,  1975  at  Taskinas  Creek;  portions 
of  this  are  shown  in  Figures  25  to  30.  The  figures  depict  incoming 
solar  radiation  in  Taskinas  Creek,  tidal  fluctuation,  and  pH  and  tem- 
perature at  the  surface  and  10  cm  below  the  sediment  surface  in  several 
marsh  zones.  These  figures  readily  show  the  high  degree  of  variability 
present  in  the  marsh,  and  also  the  difficulty  in  analyzing  this  type  of 
continuous  data.  The  curves  of  solar  radiation  (Figure  25)  and  tidal 
height  (Figure  26)  were  relatively  smooth  and  easy  to  interpret,  since 
the  diurnal  periodicity  of  radiation  and  semidiurnal  nature  of  the  tides 
were  readily  apparent.  However,  graphical  analyses  of  temperature  and 
pH  curves  is  a formidable  task,  since  patterns  of  behavior  are  masked 
by  wide  variation  and  signal  noise  (Figures  27  to  30) . 


AD-A051  929  VIRGINIA  UNIV  CHARLOTTESVILLE  DEPT  OF  ENVIRONMENTAL  —ETC  F/G  6/6 

MODEL I N6  OF  ECOLOGICAL  SUCCESSION  AND  PRODUCTION  IN  ESTUARINE  M— ETC (U) 
NOV  77  J C ZlEMANr  V E ODUM  DACV39-73-C-0121 

UNCLASSIFIED  VES-TR-D-77-35  NL 


METERS 


2000  2400 


1200 


2000  2400 


HOURS 


Figure  26.  Tidal  height  at  Taskinas  Marsh  from 
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Figure  27.  Temperature  at  the  sediment  surface  at  Taskinas  Marsh  from 
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Figure  30.  pH  at  10  cm  in  the  sediment  at  Taskinas  Marsh  from  25  to  28  Nc 


iskinas  Marsh  from  25  to  28  November  1975. 


102.  Because  the  data  were  continuous  and  spanned  several  cycles 
of  the  driving  function,  they  could  be  analyzed  by  rigorous  mathemati- 
cal techniques,  including  Fourier  series  analysis  (Rayner,  1971;  Platt 
and  Denman,  1975).  This  technique  is  utilized  to  determine  if  an  envir- 
onmental forcing  function  such  as  tide  or  radiation  exerts  an  effect 

on  another  parameter.  Alternatively,  it  can  be  used  to  reduce  large 
amounts  of  data,  such  as  the  voluminous  time  series  records  collected 
here.  A Fourier  series  is  a composite  of  sines  and  cosines  that  repre- 
sent the  decomposition  of  a time  series.  In  this  sense  it  is  a form 
of  curve  fitting,  because  the  sum  of  a set  of  sines  and  cosines  results 
in  the  replication  of  the  original  curve  3(t),  where: 

j-1 

S(t)  = ao  + E a sin(2irf  At)  + b. cos(2irf  At) 
i=l  1 1 1 1 

S(t)  = approximation  of  original  time  series 

ao,  a^,  b^  = Fourier  coefficients 

f ^ = i*"^1  harmonic  of  fundamental  frequency 

At  = data  sample  time  interval 

The  coefficients  ao,  a^,  and  b^  are  Fourier  coefficients  and  reflect  the 
influence  of  each  sine  and  cosine  component  in  the  reconstitution  of 
the  original  information. 

103.  For  example,  if  the  original  time  series  S(t)  is  a tidal 
record  with  zero  amplitude  at  time  zero  (t^) , the  Fourier  series  repre- 
sentation of  this  curve  would  be 


S (t)  = aisin(2irfiAt) 
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where  sin(2irt/T)  represents  a frequency  of  2 cycles/day.  All  of  the 
other  frequencies  would  have  very  small  coefficients  associated  with 
them  since  tides  in  the  Virginia  area  are  nearly  sinusoidal,  with  a pri- 
mary frequency  of  2 cycles/day  (U.S.  Dept,  of  Commerce,  1974).  The 
above  information  is  often  displayed  in  graphical  format  resulting  in 
what  is  called  a power  spectrum.  This  graph  depicts  the  relationship 
between  frequency  and  power.  The  term  power  signifies  the  importance  of 
a particular  frequency  in  reconstituting  the  original  data  and  is  a 
function  of  the  Fourier  coefficients.  Thus,  for  a tidal  curve,  a power 
spectrum  would  be  composed  of  a single  spike  at  a frequency  of  2 cycles/ 
day  (Figure  31),  all  other  frequencies  having  very  little  power. 

104.  In  many  cases  the  power  spectrum  for  a time  series  is 
dwarfed  by  a spike  at  zero  frequency.  When  this  occurs  it  is  often 
difficult  to  determine  if  any  other  frequencies  occur  as  a significant 
component  of  the  original  data.  The  zero  frequency  represents  the  mean 
of  the  data  series,  and  a filtering  procedure  is  often  utilized  to  elim- 
inate such  low  frequencies.  Following  this  filtering,  it  can  be  ascer- 
tained whether  the  ecologically  important  frequencies  such  as  1 and  2 
cycles/day  are  present  in  the  time  series  under  study.  For  example, 
consider  a pH  time  series  for  surficial  sediment  in  a salt  marsh.  The 
pH  may  vary  around  a mean  of  approximately  5 to  7 pH.  To  analyze  this 
time  series  for  frequencies  in  the  range  of  1 to  2 cycles/24  hours,  the 
mean  must  be  removed  from  the  data.  This  is  often  accomplished  by  util- 
izing a differencing  filter,  or  simply  by  subtracting  the  mean  of  the 
data  set  from  each  individual  data  point.  In  the  case  of  the  differ- 


encing filter  the  data  point  of  a time  series  at  time  t is  subtracted 


i 


from  the  data  point  at  time  t+1,  thus  constituting  a new  value  at  time 
t.  The  process  is  continued  for  the  entire  time  series.  The  result  of 
this  mathematical  manipulation  is  the  loss  of  low  frequencies  from  the 
original  time  series. 

105.  The  power  spectrum  analyses  were  conducted  on  continuous  moni- 
toring data  for  which  a sufficient  number  of  data  points  existed.  Power 
spectra  for  tide,  radiation,  pH,  and  temnerature  at  the  surface  of  the 
salt  marsh  are  displayed  in  Figures  31-36.  These  data  were  collected  at 
low  marsh  and  high  marsh  sites  at  Taskinas  Creek  for  the  period  of  25  to 
28  November,  1975.  The  tide  and  radiation  data  were  unfiltered,  while  the 
temperature  and  pH  data  were  manipulated  by  the  use  of  a differencing  fil- 
ter. As  might  be  expected,  the  power  spectra  for  the  tide  and  radiation 
data  have  spikes  at  0.089  and  0.044  cycle/hour,  respectively  (two  cycles 
and  one  cycle,  respectively,  per  24  hours).  There  is  little  power  at  any 
other  frequency  except  zero  since  the  mean  was  not  subtracted  from  these 
time  series.  The  power  is  primarily  at  one  frequency  since  tide  and  radi- 
ation are  both  predominantly  single-frequency  forcing  functions. 

106.  If  the  pH  and  temperature  time  series  data  were  analyzed  to 
determine  the  effect  of  the  tide  and  radiation  forcing  functions,  the 
resultant  power  spectra  for  pH  (Figure  33)  indicate  this  parameter  in 
the  low  marsh  appeared  to  have  a significant  frequency  component  around 
0.089  cycle/hour  (tidal  frequency),  while  this  frequency  was  absent 

in  the  high  marsh  pH  data.  The  surficial  temperature  data  also 
reflect  the  tidal  influence  elucidated  above.  As  is  seen  in  Figure  34, 
the  temperature  power  spectrum  for  low  marsh  had  a large  spike  at  the 
tidal  frequency  of  0.089  cycle/hour.  The  temperature  time-series  data 
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Figure  34.  Power  spectrum  of  temperature  at  the  surface  in  the 
Spartina  alterniflora  zone  at  Taskinas  Marsh. 
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Figure  35.  Power  spectrum  of  pH  at  the  surface  in  the  Spartina  patens 


for  the  high  marsh  lacked  the  tidal  frequency  but  had  a large  frequency 
component  near  1 cycle/day,  which  was  the  primary  radiation  frequency. 

107.  This  analysis  supplements,  and  to  a degree  substantiates, 
some  of  the  earlier  conclusions  derived  from  the  nutrient  analysis  of 
marsh  sediments.  It  indicates  that  the  low  marsh  was  being  affected 
to  a greater  extent  by  tidal  forces  than  the  high  marsh. 
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PART  VII:  MARSH  SUCCESSION  MODEL 

108.  One  of  the  major  objectives  of  this  project  was  to  develop 
a model  to  simulate  the  growth  and  succession  of  major  plant  species 
on  natural  marshes  and  on  dredged  material  disposal  sites.  Accordingly, 
model  development  paralleled  the  field  phase  of  the  project.  To  pro- 
duce a simulation  model  of  plant  growth  and  succession  on  salt  marshes, 
several  modeling  approaches  may  be  taken.  The  ideal  approach  is  to 
work  from  a complete  understanding  of  the  underlying  ecological  prin- 
ciples that  cause  the  growth  and  death  of  plants.  This  yields  a model 
that  can  accurately  depict  the  response  of  a plant  to  the  environ- 
mental factors — physical,  chemical,  and  biological — that  influence  the 
study  species.  However,  the  biology  of  most  organisms  is  insuffi- 
ciently understood  to  develop  models  at  this  level  of  resolution.  The 
other  extreme  in  model  design  is  a model  whose  linkages  are  totally 
empirical.  Such  a device  is  usually  developed  to  test  the  accuracy  of 
predicting  outputs  from  a series  of  chosen  inputs,  with  little  atten- 
tion given  to  the  realism  of  the  model's  internal  processes.  Super- 
ficially, this  type  of  model  seems  ideal  for  many  purposes.  However, 
as  any  type  of  model  is  produced,  it  is  inevitably  necessary  to  make 
some  simplifying  assumptions  ; and  in  the  case  of  the  highly  empirical 
model,  the  assumptions  may  severely  limit  the  usefulness  of  the  final 
product . 

1C9.  Early  experimental  models  produced  for  this  study  were 
highly  realistic  growth  models  based  primarily  on  studies  involving 
Spartina  alternlf lora.  However,  their  extreme  complexity  made  them 
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difficult  to  use.  One  of  the  major  concerns  in  this  project  was  that 
the  final  model  should  be  easily  manlpulable  for  individuals  with  only 
a basic  understanding  of  the  use  of  computers;  it  was  therefore  judged 
that  the  highly  realistic  models,  which  demand  thorough  acquaintance 
with  both  computer  science  and  salt  marsh  ecology,  were  unsuitable. 

110.  The  inherent  complexity  of  the  initial  models  also  empha- 
sized another  problem  that  aided  in  the  development  of  the  present 
model.  Like  any  complex  mechanical  device,  a model  is  limited  by  its 
weakest  components  (consequently  constructing  a model  that  is  complex 
in  some  of  its  functions  but  crude  in  others  can  often  be  futile  exer- 
cise). The  marsh  literature,  though  extensive,  is  largely  descriptive 
in  nature,  and  little  information  is  sufficiently  quantitative  to  be 
applied  toward  development  of  a simulation  model.  Spartina  alterni- 

f lora  has  been  studied  extensively,  yet  there  are  sufficiently  large 
uncertainties  about  the  physiological  ecology  of  this  plant  that  con- 
structing a realistic  growth  model  is  difficult.  Knowledge  of  the 
physiology  and  ecology  of  S.  patens  and  _D.  spicata  is  even  more  limited. 
Thus,  in  constructing  a realistic  model  using  these  three  species,  one 
must  make  a choice:  either  be  constrained  by  the  species  for  which 
the  least  information  is  available,  or  make  assumptions  about  this 
species'  responses  based  on  the  behavior  of  closely  related  species. 

111.  Since  assumptions  must  always  be  made  in  building  a model, 
constraints  must  be  imposed  on  the  model's  use.  For  instance,  if  a 
growth  model  is  developed  to  simulate  the  response  of  a species  to  sev- 
eral factors  within  an  assumed  salinity  range  of  10  to  35  percent,  a 
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salinity  input  of  5 percent  may  produce  a response  that  is  seriously 
misleading. 


Model  Concept 


112.  The  central  units  of  the  succession  simulation  were  sub- 

* 

models  of  the  growth  dynamics  of  each  species,  called  species  modules. 
The  original  species  module  developed  was  for  Spartina  alternif lora, 
because  there  was  more  information  available  for  this  species  than  for 
other  salt  marsh  plants.  A frequent  drawback  of  simulation  models 
developed  with  scarce  data  was  their  inability  to  evolve  as  the  data 
base  increased.  The  species  modules  were  developed  to  be  open-ended; 
that  is,  as  more  information  became  available  documenting  the  response 
of  the  plants  to  environmental  factors,  it  could  be  used  to  increase 
the  realism  of  the  module. 

113.  The  same  basic  species  module  structure,  considered  to  be 
representative  of  a generalized  salt  marsh  plant,  was  used  for  all  spe- 
cies (Figure  37).  The  differing  responses  of  the  individual  species 
were  incorporated  by  imbedding  response  curves  within  each  module  to 
characterize  the  response  of  given  species  to  selected  environmental 
variables. 

114.  The  model  described  here  utilized  three  species  modules, 
representing  j3.  alternif lora,  _S.  patens,  and  spicata,  but  possessed 
the  capability  of  utilizing  up  to  10  species  modules  whether  they  were 
used  to  represent  individual  species  or  functional  groups  of  species. 


115.  The  species  modules  were  combined  into  a point  model  that 


Figure  37.  Conceptual  diagram  of  a single  species  module.  The  curves  represent  a generalized 
response  of  the  plant  (y-axis)  and  the  level  of  the  driving  variable  (x-axis) . 


simulates  the  growth  of  each  species  on  a particular  square  meter  of 
marsh  surface  (Figure  38).  Because  each  point  model  contains  all  species 
modules,  it  can  simulate  the  growth  of  co-occurring  species  at  a parti- 
cular point  in  the  marsh  through  time.  The  current  modal  can  depict 
plant  growth  and  succession  on  a transect  across  a salt  marsh  by  using 
physical  variables  to  link  series  of  point  models.  Coupling  of  the 
point  models,  each  of  which  represents  a discrete  point  on  the  marsh  sur- 
face, enables  simulation  of  plant  growth  bo*-h  temporally  and  spatially 
(Figure  39)  on  the  marsh  from  the  tidal  creek  to  the  high  marsh. 

116.  The  current  version  of  the  model,  MARSH1,  contains  the  three 
species  modules  previously  described,  and  is  driven  by  four  environ- 
mental factors:  solar  radiation,  temperature,  soil  salinity,  and  tidal 
inundation.  A wide  variety  of  mathematical  formulations  are  available 
to  describe  the  growth  of  the  plant  using  such  a set  of  driving  vari- 
ables. Probably  the  most  common  method  utilizes  a multiplication  tech- 
nique where  growth  is  a result  of  some  factor  "a"  times  a factor  "B" 
times  a factor  "C"  and  so  on  so  that  the  result  is  a product  of  a series 
of  driving  variables.  This  is  often  a very  unrealistic  way  to  simulate 
a response  to  a series  of  variables  as  it  assumes  interactions  among 
the  variables  that  may  be  ill-understood  or  nonexistent. 

117.  Plant  physiologists  have  clearly  established  the  interac- 
tion of  temperature  and  solar  radiation  in  controlling  photosynthetic 
rate.  The  curves  of  Gates  (1968)  depict  photosynthesis  as  a function 
of  various  light  levels  and  differing  temperatures.  However,  detailed 
information  is  not  available  for  most  marsh  plants;  thus  it  was  felt 
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ire  38.  Conceptual  diagram  of  the  linking  of  several  species  modules  with  environmental 
forcing  functions  to  form  a model  of  a single  point  in  the  marsh. 
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Figure  39.  Conceptual  diagram  showing  the  coupling  of  several  point  models 
to  produce  a coupled  model  of  transect  with  several  stations. 
The  middle  and  lower  sections  represent  generalized  plant 
species  distributions  on  the  marsh.  The  dashed  line  in  the 
lower  section  represents  mean  high  water. 


that  a better  approach  would  be  to  devise  a series  of  independent  limit- 
ing functions  formulated  according  to  Liebig's  concept  that  implies 
that  of  the  variety  of  factors  which  control  the  growth  of  organisms, 
there  is  one  element  in  least  supply  at  any  given  time,  termed  the 
"limiting  factor." 

118.  At  each  time  step  the  model  calculates  the  amount  of  growth 
that  is  possible  due  to  each  of  the  independent  driving  parameters. 

The  model  then  checks  to  determine  which  of  these  calculated  functions 
is  the  most  limiting,  or  conversely,  which  permits  the  least  growth. 

The  incremental  growth  function  is  then  determined  according  to  the 
value  of  this  limiting  parameter.  It  should  also  be  noted  that  the 
model's  driving  variables,  which  will  be  described  in  detail  below,  are 
primarily  physical  in  nature.  This  is  due  in  part  to  the  paucity  of 
information  currently  available  on  direct  effects  of  biological  factors, 
such  as  competition,  upon  marsh  plant  growth.  The  model  has  been 
designed  so  that  other  factors  that  influence  plant  growth,  such  as 
soil  nutrients,  may  be  added  as  sufficient  information  becomes  available . 
Interactions  between  factors,  both  positive  and  negative,  may  also 
be  added  when  these  are  sufficiently  understood. 

t 

Model  Structure 

119.  Using  these  concepts,  the  model  MARSH1  was  formulated,  and 
its  general  structure  is  shown  in  Figure  40.  When  a set  of  data  is  read 
into  the  program,  it  is  immediately  printed  out  with  appropriate 
descriptions  and  formatting.  The  inputting  data  fullv  fall  into  three 


general  categories:  site-specific  data,  model  parameterization  data, 
and  plant  parameterization  data.  The  site-specific  data  include  the 
numbers  of  stations  along  the  transect,  the  number  and  species  of 
plants  to  be  included  in  the  parameterization,  and  the  latitude  of  the 
site.  The  model  parameterization  values  include  such  information  as 
the  type  of  growth  algorithm,  the  time  step  for  the  numerical  integra- 
tion, and  intervals  for  printing  and  plotting.  The  plant  parameteriza- 
tion section  includes  plant-specific  intrinsic  growth  rates  and  the 
curves  of  species'  response  to  the  driving  variables.  The  most  fre- 
quently changed  cards,  the  site-specific  data,  are  placed  first  in  the 
data  deck,  in  order  to  facilitate  manipulation. 

120.  Next,  the  time  increment  is  set  and  the  subroutine  SSPECIE 
is  called.  A simplified  flow  chart  of  SSPECIE  is  included  in  Figure 
40  but  will  be  explained  in  more  detail  in  the  following  paragraphs. 
When  all  plant  growth  values  for  a particular  computer  run  are  calcu- 
lated the  program  returns  to  the  main  flow  and  the  results  of  the  simu- 
lation are  printed  or  plotted  with  an  internal  plot  routine,  or  both. 

121.  The  bulk  of  the  computation  is  performed  within  the  sub- 
routine SSPECIE.  The  flow  diagram  for  this  subroutine  is  shown  in 
Figure  41, 

122.  Upon  entering  the  subroutine,  the  initial  procedures  set 
the  counters  of  two  loops  to  proper  initial  values.  The  general  flow 
is  through  the  two  DO  loops — the  inner  loop  (KO)  controls  the  parti- 
cular species  of  plant  being  dealt  with  at  the  moment,  and  the  outer 
loop  (NN)  determines  the  particular  station  on  the  marsh  for  which 
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Figure  41.  Flow  diagram  of  the  main  subroutine,  SSPECIE. 


plant  biomass  Is  being  calculated.  The  routine  begins  at  station  1 
(NN=1)  and  at  this  station  calculates  the  biomass  of  each  species  (K)-L 
to  NUM  (where  NUM  = the  n**1  species  under  consideration) , based  on  the 
minimum  limiting  coefficient. 

123.  The  first  function  in  the  subroutine  SSPEC1E,  determines  the 
level  of  solar  radiation  incident  upon  the  plant  species  under  consid- 
eration. The  basis  for  the  driving  curves  (Figure  42)  were  seasonal 
values  calculated  from  twenty-year  means  given  in  the  Climatological 
Atlas  (U.S.  Dept,  of  Commerce,  1973).  These  averaged  seasonal  values 
for  light  along  the  the  eastern  seaboard  between  25°  to  45°N  show 
graphically  how  light  may  act  as  a control  in  plant  growth  along  a lati- 
tudinal gradient.  The  maximum  values  reached  during  the  summer  are 
quite  consistent  throughout  the  entire  range  shown,  suggesting  that 
light  would  not  be  any  more  limiting  in  the  higher  latitudes  in  the 
summer  than  in  the  lower  latitudes.  However,  the  minimum  winter  light 
levels  are  dramatically  different.  At  30°N,  roughly  the  latitude  of 
Jacksonville,  Florida,  the  light  levels  in  December  are  as  high  as  the 
level  at  40°N,  near  Delaware  Bay  and  the  New  Jersey  marshes,  in  early 
October  (vigure  42). 

124.  The  curves  shown  are  incorporated  in  the  model  as  a set  of 
third-order  polynomial  equations.  The  output  values  represent  sunlight 
intensity  at  the  top  of  the  plant  canopy.  If  a selected  latitude  is 
intermediate  between  any  two  of  the  five  curves,  a routine  calculates 
an  interpolated  value. 

125.  Competition  between  plant  species  at  a site  is  determined 
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by  FUNCTION  COHT,  which  assumes  that  all  individuals  of  a particular 
species  at  a site  are  equal  in  height.  The  height  differentials  among 
species  are  utilized  to  calculate  the  ambient  light  level  among  the 
plants  by  simulating  attenuation  of  sunlight  from  the  canopy  to  the 
sediment  surface.  The  amount  of  radiation  as  a function  of  height 
within  the  stand  is  calculated  in  Function  COSHADE  (see  Figure  41) . 

It  Is  important  to  note  that  due  to  differences  in  the  density  of  foli- 
age, different  plant  zones  in  the  marsh  have  varying  degrees  of  radia- 
tion attenuation  with  depth  of  canopy  (Gates,  1968). 

126.  The  next  environmental  factor  considered  is  temperature. 

Again  the  Climatological  Atlas  (U.S.  Dept,  of  Commerce,  1973)  was  used 
to  obtain  data  on  the  seasonal  variation  of  temperature  as  a function 
of  latitude.  At  the  lower  latitudes,  temperatures  show  little  seasonal 
change  up  to  27°N,  varying  from  15-20°C.  In  contrast,  seasonal  tempera- 
tures indicate  that  vegetation  in  the  high  latitudes  experiences  a wide 
seasonal  range  of  temperature  with  a lower  mean  yearly  maximum  (16-22°C) 

(see  Figure  43).  As  with  the  radiation  data,  this  information  is  incor- 
porated into  the  model  with  third-order  polynomial  equations.  These 
formulations  are  embedded  in  Function  TEMP  (see  Figure  41).  When  lati- 
tude data  are  entered,  the  monthly  variation  of  temperature  at  the  loca- 
tion specified  will  be  determined. 

127.  One  of  the  most  important  environmental  parameters  affecting 
plant  zonation  is  tidal  inundation.  This  has  been  recognized  by  numer- 
ous investigators  and  was  further  substantiated  in  other  sections  of 

- 

this  report.  It  was  shown  herein  that  tidal  inundation  is  the  most 
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Figure  43.  Seasonal  temperature  values  for  the 
eastern  coast  of  the  United  States. 


significant  factor  affecting  the  levels  of  marsh  soil  nutrients  and  the 
resultant  marsh  plant  productivity. 


128.  A computer  program.  Program  TIDE,  was  developed  at  the  Uni- 
versity of  Virginia  to  calculate  the  time  of  tidal  inundation  and  the 
number  of  tidal  inundations  per  year  at  any  height  within  a marsh.  The 
data  base  for  this  program  includes  times  and  heights  of  each  daily 
high  and  low  tide  at  a series  of  stations  along  the  East  Coast  from  Key 
West  to  Maine.  The  results  of  this  analysis  (Figure  44)  indicate  that, 
with  ninor  variations,  each  of  the  curves  can  be  resolved  into  two 
nearly  linear  slopes.  The  lowest  curve  on  Figure  44,  the  inundation 
curve  for  the  Patuxent  River,  serves  for  illustration.  The  highest 
point  inundated  by  maximal  spring  tides  is  about  1.01  m above  mlw.  The 
area  from  mean  high  tide,  approximately  0.76  m,  to  this  point  receives 
minimal  inundation,  as  shown  by  the  slope  of  the  curve  at  this  point. 
However,  the  inundation  time  between  mlw  and  mhw  is  an  essentially  lin- 
ear function  from  0 to  0.76  m.  Thus,  while  the  total  tidal  range  from 
mlw  to  spring  high  water  is  1.01  m,  the  upper  25  percent  of  the 
tidal  range  receives  less  than  5 percent  of  maxir  lm  inundation 

129.  While  all  of  the  inundation  curves  depicted  show  a similar 
form,  the  pattern  is  extended  in  regions  of  extremely  large  tidal  ranges^ 
such  as  Maine  and  parts  of  Georgia.  The  relationship  between  the  maxi- 
mum tide  and  the  inflection  point  is  shown  in  Figure  45  and  indicates 
the  consistency  of  the  relationship.  Because  of  the  high  ;orrelation 
between  the  inflection  point  and  the  maximum  tidal  height,  the  series 

of  curves  for  tidal  inundation  can  be  condensed  into  a normalized 
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Figure  45.  Relationship  between  peak  tidal  height  (maximum  predicted 
spring  tide)  and  the  projected  point  where  the  main  body 
of  the  curve  would  intercept  the  vertical  axis. 
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relationship  (Figure  46)  by  taking  each  curve  from  Figure  44  separately 
and  dividing  the  value  on  each  axis  by  the  maximum  value  for  that  site 
relating  both  the  height  and  inundation  values  of  percent  of  minimum. 

The  normalized  curve  (Figure  46)  permits  calculation  of  the  time  of 
inundation  given  the  height  of  the  station  of  interest  and  height  of 
the  spring  tides  above  mlw.  The  inflection  point  in  the  normalized 
curve  occurs  at  75  percent  of  the  maximum  tidal  height  and  corresponds 
to  mhw.  This  is  also  the  level  on  a marsh  where  the  interface  between 
the  j5.  alternif lora  zone  and  j3.  patens  - D.  spicata  zone  usually  occurs. 
The  region  above  mhw  is  inundated  less  than  4 percent  of  the  year. 

130.  Within  the  succession  model,  the  relationship  between  rela- 
tive height  and  percent  inundation  is  represented  by  a piecewise  linear 
formulation  embedded  in  subroutine  TIDE. 

131.  The  final  environmental  parameter  considered  is  soil  salin- 
ity. Figure  12  shows  changes  in  soil  salinity  along  a series  of  tran- 
sects in  Bennett's  and  Taskinas  marshes.  Salinity  in  the  hydrologically 
active  low  marsh  is  comparable  to  that  found  in  the  contiguous  marsh 
creek.  Higher  salinity  is  found  at  higher  elevations  above  mlw.  This 
trend  is  simulated  in  the  model  by  a piecewise  linear  formulation 
(Figure  47),  which  is  represented  by  graphing  relative  height  in  a 
marsh  (station  height/mean  tidal  range)  versus  salinity.  This  graph 
shows  that  an  increase  in  elevations  on  the  marsh  surface  results  in  a 
soil  salinity  nearly  double  that  found  in  the  nearby  marsh  creek;  the 
difference  may  be  caused  by  evaporation  and  possibly  by  water  loss  due 
to  transpiration.  Where  the  marsh  borders  the  upland,  freshwater  run- 
off dilutes  soil  salinity  to  almost  0 ppt. 
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132.  In  order  to  determine  how  much  each  species  is  curtailed  by 


environmental  interaction,  a series  of  parameterization  curves  were 
developed.  These  curves,  mathematically  embedded  within  the  model, 
represent  plant  response  to  solar  radiation,  temperature,  tidal  inunda- 
tion, and  soil  salinity.  The  x-axis  of  the  parameterization  curve 
reflects  the  level  of  the  chosen  forcing  function,  while  the  y-axis 
represents  a normalized  growth  response.  The  response  of  a plant  to  a 
specific  level  of  an  environmental  forcing  function  ranges  between  0 
and  1;  if  the  plant  response  is  0,  the  environmental  factor  has  inhi- 
bited growth;  a response  of  1 indicates  the  level  of  the  factor  is  opti- 
mum  for  growth. 

133.  Figure  48  depicts  the  responses  of  S.  alterniflora  to  temp- 
erature, radiation,  percent  inundation,  and  salinity.  The  radiation  and 
temperature  parameterization  curves  are  similar  to  those  presented  by 
Gates  (1968) . These  curves  represent  generalized  plant  responses  to 
the  above  parameters,  since  they  were  generated  from  data  collected  on 

» many  phylogenetically  divergent  species.  The  response  of  JS.  alterni- 

flora to  salinity  reflects  the  fact  that  this  plant  is  considered  a 
facultative  halophyte.  Therefore,  an  increase  in  salinity  results  in  a 
decreased  biomass  response.  The  percent  inundation  parameterization 
curve  indicates  the  tolerance  of  this  plant  to  submergence  since  this 
plant  is  found  only  in  the  low  marsh. 

134.  The  response  of  j>.  patens  and  D.  spicata  to  each  of  the 

environmental  parameters  is  depicted  in  Figure  49.  Competition  between 

■ i 

the  species  is  not  modeled  explicitly;  however,  the  data  do  reflect 

natural  situations.  Lacking  data  to  the  contrary,  the  parameterization 
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Figure  48.  Normalized  growth  response  of  Spartina  alterniflora  to 
environmental  driving  parameters. 
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Figure  49.  Normalized  growth  response  of  Spartina  patens  and  Dis- 
tichlis  spicata  to  environmental  driving  parameters. 
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curves  for  radiation  and  temperature  are  assumed  to  be  identical  to 
those  of  j>.  alternif lora.  However,  the  actual  responses  of  J5.  alterni- 
flora  and  S.  patens  are  not  identical,  since  their  heights  and  the 
resultant  shading  effects  are  usually  quite  different.  It  was  assumed 
that  J5.  patens  is  better  adapted  to  higher  salinities  since  it  is  found 
in  the  high  marsh.  This  is  reflected  in  the  salinity  parameterization 
curve.  Similarly,  as  a result  of  its  position  in  the  marsh,  maximum 
growth  response  of  j>.  patens  to  percent  inundation  was  found  at  a lower 
value  than  that  of  S^.  alternif  lora. 

135.  D.  spicata  and  j>.  patens  have  similar  responses  to  the  four 
environmental  parameters  (Figure  49).  The  difference  lies  in  the  fact 
that  the  tolerance  to  salinity  and  percent  inundation  of  D.  spicata  is 
greater  than  that  of  S^.  patens.  Parameterization  curves  for  j>.  alterni- 
f lora,  S_.  patens,  and  D.  spicata  have  been  discussed  because  these 
three  plants  dominate  the  marsh  vegetation  in  many  localities,  but 
curves  for  any  other  species  can  be  formulated  if  the  necessary  field 
data  are  available. 

136.  At  this  point  in  the  program  a plant  response  value  or 

1 

limiting  coefficient  is  determined  for  each  species  and  each  environ- 
mental factor.  The  assumption  that  certain  environmental  parameters 
are  largely  responsible  for  restricting  growth  (law  of  the  minimum) 
simplifies  the  task  of  determining  the  relative  effect  of  the  limiting 
factors.  At  each  calculation  the  limiting  factor  is  stored,  so  that 
the  operator  may  recall  at  any  time  the  parameter  that  is  limiting  a 
particular  species.  With  this  information,  a Runge-Kutta  numerical 
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I integrating  subroutine  is  accessed  to  calculate  plant  productivity. 

Within  this  subroutine  two  options  for  calculating  productivity  are 
available:  a logistic  method  and  a maximum  productivity  method.  Several 
assumptions  are  inherent  in  using  these  methods.  In  the  logistic  for- 
mulations productivity  is  determined  as  a function  of  the  intrinsic 
growth  rate,  the  limiting  coefficient,  and  the  contemporaneous  biomass. 
Also  included  is  a negative  feedback  term  associated  with  environmental 
inhibition.  The  maximum  productivity  algorithm  is  the  product  of  the 
limit  .ng  coefficient  multiplied  by  maximum  productivity  per  unit  time. 
Within  this  model  alterations  are  made  for  relative  elevation  and  all 
species  of  interest.  The  output  yields  plotted  and  printed  biomass 
estimates  for  each  species  as  a function  of  time  and  space. 

137.  As  described  above,  the  marsh  model  is  composed  of  a series 
of  point  models  that  are  linked  spatially  (Figure  40).  The  real  world 
basis  of  this  linkage  is  the  lateral  transport  of  water-borne  sediments. 
Within  the  model,  relative  elevation  in  the  marsh  controls,  to  a great 
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extent,  the  interaction  between  plant  and  environment.  At  each  loca- 
tion in  the  marsh,  allowance  is  made  for  a temporal  increment  of  sedi- 
ment, and  the  resulting  changes  in  elevation  control  the  spatial  dis- 
tribution of  plant  species. 

138.  The  above  discussion  serves  to  illustrate  the  logic  and 
assumptions  inherent  in  the  succession  model.  A complete  listirg  of 
the  model,  with  detailed  instructions  on  assembly  of  the  data  cards  and 
use,  is  contained  in  Appendix  B.  Using  only  the  parameterization 
curves  given  for  the  three  species  utilized,  a vast  number  of  simula- 
tions are  possible.  The  following  discussion  will  illustrate  the 
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results  of  one  set  of  simulations. 


139.  For  these  experimental  simulations,  latitude  and  tidal 
range  were  varied  to  study  the  effect  of  elevation  on  plant  growth  and 
species  location  on  the  marsh  surface.  The  first  simulation  (Figures 
50-55),  the  growth  of  plants  on  salt  marsh  at  latitude  30°N,  with  a tidal 
creek  salinity  of  20  o/oo,  a spring  tidal  range  of  1.83  m,  and  a mean 
tidal  range  of  1.3  m was  depicted  at  stations  varying  from  0.31  to  1.33  m 
above  mlw.  The  three  letters  at  the  right-hand  side  of  the  figure  illus- 
trate which  of  the  four  environmental  variables  was  limiting  for  each 
species  (T  = temperature,  H = tidal  height,  S = salinity,  and  L = light). 
The  first  letter  refers  to  the  environmental  variable  limiting  the  first 
species  listed,  the  second  letter  refers  to  the  second  species,  and  the 
third  letter  refers  to  the  third  species. 

14C.  _S.  alternif lora  was  found  growing  in  the  low  marsh  as  ex- 

pected with  growth  beginning  at  0.62  m above  mlw  or  near  the  mean  tidal 
height.  With  an  increase  in  elevation,  the  pure  stand  of  j>.  alternif lora 
graded  into  a mixed  association  including  S_.  patens  and  D.  spicata  near 
the  mean  high  tide  mark.  On  Figures  50-55,  the  graphs  indicate  species 
1 (S.  alterniflora)  and  3 (D.  spicata)  only.  S^.  patens  (#2)  was  present 
but  had  the  same  value  as  I).  spicata.  When  this  occurs,  the  plotting 
routine  prints  the  single  highest  species  number  encountered.  In  the 
high  marsh  above  mlw  only  S^.  patens  and  D.  spicata  were  found.  In  this 
simulation,  S..  patens  and  D.  spicata  were  excluded  from  the  low  marsh  due 
to  excessive  inundation,  while  growth  on  the  high  marsh  for  these  plants 
was  most  limited  due  to  light  and  temperature. 

141.  In  the  second  simulation  shown  here  (Figures  56-60),  the 


139 


1.  non 


n*«i  Motuss 


1M4.I  MINI 
111 


Figure  50.  Results  of  a simulation  at  latitude  30  N,  with  a creek  salinity  of  20  o/oo,  a me, 
tide  range  of  4 ft  (1.22  m) , a maximum  range  of  6 ft  (1.83  m) , and  a time  period 
of  48  months.  This  first  figure  depicts  the  simulation  at  the  1-ft  (0.31-m)  ele 
vation  above  ralw  on  the  marsh. 


Figure  52.  Depicts  the  simulation  at  the  3-ft  (0.91-m)  elevation  above  mlw  on  the  marsh. 


gure  55.  Depicts  the  simulation  at  the  6-ft  (1.83-m)  elevation  above  mlw  on  the  marsh 


Figure  56.  Results  of  a simulation  at  latitude  35  N,  with  a creek  salinity  of  20  o/oo, 
a mean  tide  range  of  4 ft  (1.22  m) , and  a maximum  range  of  6 ft  (1.83  m) 
for  48  months.  This  figure  depicts  the  simulation  at  the  2-ft  (0.61-m) 
elevation  above  mlw  on  the  marsh. 
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Figure  58.  Depicts  the  simulation  at  the  4-ft  (1.22-m)  elevation  above  mlw  on  the  marsh 
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marsh  was  located  at  a latitude  of  35  N,  and  all  other  parameters  were 
held  constant.  The  location  of  species  within  the  marsh  with  respect 
to  tidal  height  was  similar  to  that  for  the  marsh  at  30°N,  but  the  tem- 
poral change  of  biomass  and  peak  biomass  varied  considerably.  In  this 
more  northerly  area,  the  growing  season  would  be  shorter.  This  was 
reflected  in  a more  rapid  increase  in  biomass  after  growth  is  initiated. 
A slightly  lower  peak  biomass  level  was  found  at  this  higher  latitude, 
but  because  the  growing  season  is  shorter,  the  total  area  under  the 
curve  is  considerably  less  than  at  30°N.  This  relationship  was  found 
for  all  locations  within  the  marsh. 

142.  The  simulations  used  for  illustration  showed  that  radia- 
tion, temperature,  and,  particularly,  tidal  inundation  controlled  plant 
growth  and  plant  location  on  the  marsh  surface.  Salinity  was  seldom 
seen  to  be  the  limiting  factor  which  generated  the  above  results 
(Figures  50-60).  This,  however,  may  be  largely  due  to  an  artifact  of 
the  parameterization  of  the  model. 

v 143.  No  model  is  any  better  than  the  data  and  constructional 

information  used  to  generate  the  model.  The  model  presented  is  of  more 
utility  than  the  results  presented  here.  The  results  are  limited  by 
the  parameterizations  that  were  available  when  this  model  was  origin- 
ally constructed.  However,  the  model  was  not  so  structured,  and  is 
limited  only  by  the  knowledge  available  when  it  is  utilized.  Therefore, 
all  new  data  and  theory  generated  that  affects  growth,  distribution, 
and  succession  on  a marsh  may  be  incorporated  in  future  simulations. 
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PART  VIII:  SUMMARY  AND  CONCLUSIONS 


144.  At  the  outset  of  this  project  three  parallel  studies  were 
designed  to  gain  insight  into  the  processes  that  determine  plant  growth 
and  succession  on  estuarine  salt  marshes.  The  first  was  concerned  with 
both  the  daily  and  seasonal  ranges  of  marsh  physicochemical  parameters 
and  the  subsequent  biotic  responses  within  different  vegetative  zones. 

A series  of  marshes  were  visited  and  surveyed,  but  effort  was  concentra- 
ted on  two  intensive  sites  in  Virginia. 

145.  The  second  study  involved  obtaining  continuous  records  of 
several  selected  physical  and  chemical  parameters  from  within  the  vege- 
tation zone.  This  required  the  development  of  a continuous  recording 
data  monitoring  system  suitable  for  use  in  the  estuarine  system. 

146.  The  third  study  of  the  project  was  the  development  of  a sim- 
ulation model  of  plant  growth  and  succession  on  the  marsh.  The  initial 
models  were  developed  using  information  available  in  the  literature  on 
the  functions  of  estuarine  salt  marshes.  As  the  project  progressed, 
information  was  added  from  field  studies  and  contact  with  other  scien- 
tists in  the  field. 

147.  Data  from  the  first  study  were  analyzed  with  correlation 
analysis  of  plant  growth  related  to  habitat  condition  and  showed 
clearly  that  the  growth  of  J5.  alternif lora  exhibited  a strong  positive 
correspondence  with  elevation  and  tidal  inundation  at  both  major  study 
sites.  In  addition,  this  plant  demonstrated  a negative  correlation 
with  soil  salinity,  while  the  standing  crop  was  significantly  corre- 
lated with  iron  availability. 
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148.  For  j5.  patens  and  D.  spicata,  the  measured  responses  to 
tides  differed  conspicuously.  At  one  site  where  the  high  marsh  was 
rarely  flooded,  the  growth  of  neither  species  was  closely  related  to 
inundation  or  elevation.  At  the  other  intensively  studied  site,  which 
receives  flooding  during  spring  tides,  there  was  a strong  correlation 
of  tide  with  j>.  patens  (positive)  and  D.  spicata  (negative)  growth. 
Wherever  regular  flooding  takes  place  the  growth  of  j>.  alternif lora  and 
j>.  patens  demonstrated  a negative  correlation  with  soil  salinity  while 
the  more  salt-tolerant  D.  spicata  was  positively  correlated  with  salin- 
ity. 

149.  The  fact  that  habitat  factors  are  closely  associated  with 
both  tide  and  plant  growth  wherever  fairly  regular  inundation  occurs 
suggests  that  inundation  acts  as  a master  variable  influencing  plant 
growth  both  directly  and  by  controlling  other  features  of  the  marsh 
habitat. 

150.  Conclusions  drawn  from  the  second  study  and  the  records 
of  continuous  field  monitors  tend  to  reinforce  the  concept  of  tidal 
Inundation  as  the  master  variable  for  explaining  plant  distribution, 
especially  in  the  region  between  the  low  tide  marsh  and  the  upper 
extent  of  spring  high  tides.  The  power  spectrum  analysis  of  continuous 
monitoring  data  from  the  marsh  with  an  infrequently  flooded  high  marsh 
provides  a good  example.  The  surface  temperatures  from  the  j3.  alterni- 
f lora  zones  showed  the  major  power  peak  to  be  at  0.089  cycle/hour, 
which  is  the  tidal  frequency  for  this  area,  while  in  the  higher  jS. 
patens  - D.  spicata  zone,  the  frequency  with  the  largest  power 


153 


component  is  that  of  nearly  1 cycle/day,  corresponding  to  the  input  of 
solar  radiation.  The  analyses  of  pH  data  further  showed  a large  spike 
of  the  tidal  frequency  in  the  low  S^.  alternif lora  marsh,  but  not  in  the 
high  marsh. 

151.  The  third  objective  of  this  project,  the  development  of  a 
marsh  successional  model,  was  designed  around  the  concept  of  species 
modules  that  depict  the  growth  dynamics  of  j>.  alternif  lora,  j>.  patens, 
and  I).  spicata. 

152.  The  species  modules  were  combined  into  a point  model  that 
simulates  the  growth  of  each  species  on  a particular  square  meter  of 
marsh  surface.  Coupling  of  a series  of  point  models  enables  simulation 
of  plant  growth  both  temporally  and  spatially  across  the  marsh.  The 
total  marsh  model  is  driven  by  four  environmental  variables:  solar 


radiation,  temperature,  soil  salinity,  and  tidal  inundation.  A special 
subprogram  was  developed  to  calculate  the  time  of  tidal  inundation  and 
number  of  tidal  inundations  per  year  at  any  height  within  a marsh. 

153.  Experimental  simulations  demonstrated  that  radiation,  tem- 
perature, and,  particularly,  tidal  inundation  were  responsible  for  con- 
trolling plant  location  and  growth  on  the  marsh  surface.  Salinity,  at 
least  in  these  preliminary  runs,  seldom  appeared  to  be  a limiting 
factor.  In  the  final  analysis,  the  model  appears  to  have  great  utility 


for  future  research,  but  is  limited  at  the  moment  by  incomplete  para- 
meterizations. 
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1.  One  severe  limitation  on  the  application  of  modern  analytical 
techniques  in  estuarine  ecology  has  been  the  inability  to  measure  envi- 
ronmental variability  in  the  field  at  a rapid  rate.  Most  ecological 
measurements  are  made  at  infrequent  intervals,  and  the  time  period 
between  samples  is  relatively  long.  The  experimental  field  data  system 
described  here  was  designed  to  monitor  various  physical  and  chemical 
parameters  at  a rapid  rate  and  to  monitor  these  parameters  at  several 
sites  on  the  marsh  simultaneously. 

2.  Because  there  was  no  commercially  available  system  that  would 
perform  all  of  the  functions  desired,  a unit  was  designed  in  a modulat1- 
ized  fashion,  combining  commercially  available  components  with  sections 
whose  designs  were  tailored  to  the  project's  needs.  The  objective  of 
this  segment  of  the  study  was  to  develop  the  capability  to  compare  both 
long-term  and  short-term  variations  among  different  zones  on  the  marsh. 
For  this  reason  a system  was  designed  and  built  that  was  able  to  moni- 
tor stations  spread  across  several  hundred  meters  of  marsh  surface.  A 
similar  data  recording  system  was  constructed  by  several  other  members 
of  the  University  of  Virginia  Department  of  Environmental  Sciences  for 
T-nltoring  eutrophication  in  rivers.  The  final  report  of  this  study 
<*r' ; v et  al.,  1975)  describes  the  theory  and  circuitry  in  great  detail. 
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Configuration 

3.  Figure  A-l  illustrates  the  overall  configuration  of  the  sys- 
tem as  deployed  in  the  field,  and  Figure  A-2  depicts  the  contents  of 
one  submersible  amplifier  unit  with  sensors  and  power  supply.  The  sys- 
tem consisted  of  the  following  major  components: 

<1.  Sensors  to  measure  the  magnitude  of  the  described 
parameter. 

b^.  Cable  leads  to  connect  the  sensors  to  the  electronic 
amplifiers. 

Signal  amplifiers  with  power  sources  to  amplify  the 
measured  signal  and  connect  it  to  voltage  compatible 
with  the  data  recorders. 

<1.  Waterproof  housing  to  shield  the  amplifiers  and 
batteries. 

e.  Connecting  cables  to  transmit  the  amplified  signal  to 
the  data  recorders. 

jf.  Junction  or  distribution  boxes  to  allow  the  selection 
of  the  various  parameters  that  are  to  be  recorded  at 
a particular  time. 

£.  Digital  data  recorders  to  record  the  signals. 

h.  Waterproof  housing  to  shield  the  recorders. 


Field  Sensors 


4.  Temperature  was  measured  electronically  by  a twin  bead  ther- 
mistor and  an  accompanying  amplifier.  A Fenwal  GB  34PMG2  Iso-curve 
glass  head  thermistor  was  used.  The  variable  resistance  of  the  ther- 
mistor is  measured  by  the  amplifier  and  the  signal  converted  to  a 0.0 
to  1.0  v-d.c.  output.  The  amplifier  may  be  adjusted  so  that  the  one 
valid  output  range  corresponds  to  a variety  of  temperature  ranges.  For 
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Figure  A-l.  General  configuration  of  the  data  monitoring  system. 


Probes 


the  uses  of  this  project  the  range  selected  was  0-40°C.  The  thermistor 
is  soldered  to  a lead  wire,  and  the  thermistor  and  solder  joint  are 
molded  into  a 10-cm  length  of  clear  polymer  plastic  using  Scotchcast 
electrical  potting  compound  to  provide  electrical  insulation  and  struc- 
tural strength.  The  resulting  unit  is  exceptionally  rugged  and  resists 

damage  from  general  field  use.  Accuracy  is  approximately  ± 0.1°C. 

£H 

5.  Combination  pH/reference  electrodes  have  been  found  to  give 
satisfactory  performance  under  the  extreme  conditions  of  use  they 
would  experience  in  this  project.  Electrodes  must  withstand  submergence 
for  extended  periods  of  time  in  seawater,  and  must  be  able  to  function 
in  sediments  that  are  abrasive  and  frequently  anaerobic. 

6.  Initially  the  Thomas  4094-510  sheathed-in-glass  combination 
probe  was  used.  This  gave  satisfactory  performance,  but  was  vulnerable 
to  breakage.  Because  it  is  not  sealed,  it  required  modification  for 
this  field  use.  Furthermore,  it  has  to  be  refilled  periodically. 

7.  The  electrode  generally  used  was  the  Sensorex  Model  S300C, 

► 

which  is  a combination  electrode  with  a completely  sealed  construction. 
This  is  permanently  filled  with  a saturated  gell  KC1,  and  is  contained 
in  a high  impact  polymer  case.  These  electrodes  may  be  used  as  pur- 
chased for  short  times,  if  great  care  is  taken.  For  heavy  or  pro- 
longed use  it  proved  necessary  to  further  protect  the  cable  lead  and 
the  joint  where  the  cable  lead  enters  the  probe.  This  was  done  by 
threading  the  lead  wire  through  a length  of  1/8-in.  Tygon  tubing  and 
sealing  the  tubing  to  the  head  of  the  probe  with  household  caulking 
compound.  The  joint  was  then  securely  wrapped  with  electrical  tape. 


8.  The  pH  amplifier  was  a solid  state,  expanded  scale  pH  meter 
with  an  Analog  Devices  model  173J  Isolation  Amplifier  to  reduce  current 
leakage  problems  from  the  electrodes  and  errors  due  to  ground  loops  in 
the  highly  conductive  medium  where  the  system  must  function.  The  pH 
amplifier  gives  an  output  of  0.0  to  1.0  v-d.c.,  which  is  normally  equi- 
valent to  pH  4.0  to  10.0  ± 0.1,  although  this  range  is  adjustable. 

Eh 

9.  The  oxidation-reduction  potential  or  Eh  system  utilizes  the 
same  expanded  scale  pH  amplifier  as  the  pH  sensor,  but  as  the  signal  is 
several  orders  of  magnitude  greater  than  across  a pH  sensor  pair,  the 
isolation  amplifier  is  omitted  at  present.  The  platinum  electrode  is 
made  by  soldering  a 22  gauge  platinum  wire  about  405  cm  long  into  a 12- 
cm  brazing  rod.  This  unit  is  sealed  into  a 20-cm  pyrex  tube  with  poly- 
ester resin.  For  a reference,  this  unit  utilizes  a Sensorex  reference 
electrode  or  shares  the  reference  half  of  the  combination  pH  electrode. 
The  amplifier:  gives  0.0  to  1.0  v-d.c.  output  for  Eh  values  in  the  range 
of  ± 400  millivolts.  This  unit,  however,  has  been  the  source  of  major 
problems:  current  leakages  that  cause  erratic  Eh  recordings  and  also 
interfere  with  the  pH  readings.  Another  problem  is  that  there  are  few 
standards  with  which  to  calibrate  the  units.  The  most  commonly  used  is 

Zobell's  solution  with  an  Eh  of  +430  millivolts  at  25°C. 

Salinity 

10.  Salinity  of  the  water  in  tidal  creeks  was  measured  by  an 
Instro  Company  model  300  inductance  salinometer.  Salinity  was  deter- 
mined by  an  inductive-conductive  measurement  that  was  internally  com- 
pensated for  temperature  effects  and  linearized.  The  unit 


required  a nominal  9 to  15  v-d,c,  input  and  gave  an  output  of  0.0  to 

1.0  v-d.c.  which  corresponded  to  0 to  40  ppt. 

Solar  radiation 

11.  Incident  solar  radiation  was  measured  by  an  LI-200  SR  Pyrano- 
meter  (Lamda  Instrument  Corporation)  that  consisted  of  a silicon  solar 
cell  mounted  beneath  a diffusing  surface  and  designed  to  give  a cosine 
response.  As  the  relative  spectral  response  is  not  constant  over  the 
full  range  of  the  solar  spectrum,  the  sensor  should  be  used  only  in  the 
open  air  and  not  under  plant  canopies  or  artificial  lights.  This 
instrument  measures  radiation  directly  in  watts  per  square  meter,  giving 
total  energy  input.  Absolute  accuracy  is  ± 5 percent  when  calibrated 
against  an  Epply  Pyranometer  under  natural  daylight  and  clear  condi- 
tions. 

Tide  height 

12.  The  tidal  level  in  the  nearby  marsh  was  measured  by  a Senso- 
tec  Model  FJE  pressure  transducer,  an  absolute  pressure  unit  with  a 
reference  of  14.3  psi  sealed  in  the  reference  chamber.  It  possessed 
built-in  electronics,  received  ± 12  v-d.c.  inputs,  and  gave  an  output  of 
0.0  to  1.0  v-d.c.  corresponding  to  0.0  to  6.09  m of  water  pressure  with 
an  accuracy  of  ± 0.1  percent  of  full  scale. 


Signal  Amplifiers  and  Housings 


13.  The  sensors  for  temperature , pH , Eh,  and  solar  radiation  were 
connected  via  waterproof  cables  and  connectors  to  housing  which  con- 
tained signal  conditioning  amplifiers  for  the  sensors.  Because  of  the 


extremely  low  voltages  and  current  across  the  sensors,  especially  the 


pH  probe,  there  was  a limitation  in  the  length  of  the  lead  wire  between 
the  sensors  and  amplifiers.  In  the  field  unit,  this  length  was  kept 
below  3.5  m. 


14.  The  signal  conditioning  amplifiers  required  ± 12  v-d.c.  that 
was  supplied  by  battery  packs.  The  battery  packs  and  amplifiers  were 
housed  in  pressureproof  housings.  These  housings  were  milled  from 
10.6-cm  i.d.  PVC  tubing,  and  were  closed  with  O-ring  sealed  end  PVC 
caps.  The  battery  packs  used  sixteen  15  v-d.c.  transistor  batteries 
(D  cells)  each;  one  pack  powered  a set  of  three  or  four  amplifiers 
clustered  together  in  each  tube  for  over  two  weeks.  An  alternate 
method  to  supply  power  would  be  to  run  power  cables  from  a central  bat- 
tery source  to  each  tube  and  to  the  data  recorder.  This  option  was 
tried  initially,  but  was  discarded  due  to  problems  with  current  leaks 
and  ground  loops  between  the  sensors,  recorders,  and  amplifiers. 

15.  From  the  amplifier  tubes,  the  signals  were  fed  back  through 
waterproof  cables  and  connectors  to  a junction  box. 

16.  Certain  instruments,  such  as  the  salinometer,  current  meter, 
and  tide  gauge,  were  protected  by  their  own  self-contained  housings,  and 
linked  directly  into  the  junction  box. 

J 

Junction  Box 

17.  All  field  sensors  were  connected  to  a junction  or  distribu- 
tion  box.  This  was  currently  capable  of  receiving  the  input  from  40 
field  sensors.  A Fluke  model  8000  digital  multimeter  was  utilized  as  a 
readout  device,  and  all  of  the  40  inputs  could  be  monitored  on  the 
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digital  meter  by  dialing  the  appropriate  input  channel.  A patch  panel 
was  provided  on  the  front  of  the  junction  box,  and  any  of  the  40  input 
channels  could  be  patched  into  the  16  inputs  to  the  tape  recorder  for 
continuous  monitoring.  Thus,  16  channels  could  be  monitored  continu- 
ously by  the  recorder,  or  40  channels  of  data  might  be  monitored  inter- 
mittently by  use  of  the  digital  output  unit.  When  using  this  device 
manually,  however,  coding  was  necessary  as  the  sensor  am,  ifiers  read 
out  as  0.0  to  1.0  v-d.c.  outputs  and  not  as  scaled  values.  As  a result, 
conversion  curves  of  nomographs  were  used  to  determine  the  actual  value 
of  the  parameter  being  measured. 

Recorder 

18.  The  data  were  transmitted  from  the  junction  box  through 
moistureproof  cables  and  couplings  to  the  recorder  unit.  The  recorder 
was  a Metrodata  Dl-616  data  logger  coupled  with  a DR-99  digital  readout. 
The  data  logger  recorded  on  a 6.35-mm  endless  loop  tape  cartridge  in 
serial  format.  The  analog  input  was  converted  to  digital  format  for 
filing  on  the  tape.  Recording  times  could  be  varied  from  continuous 
(approximately  2.5  full  16  channel  scans  per  second)  to  one  scan  each 
73  minutes  via  a switch  from  a quartz  clock.  The  unit  was  capable  of 
operating  from  either  a 115  Hz  or  12  v-d.c.  source. 

19.  Although  these  units  were  rather  sophisticated,  the  field 
record  was  quite  acceptable.  The  only  consistent  stoppages  occurred 
when  the  temperature  dropped  below  freezing.  Modifications  were  re- 
quired to  decrease  the  sensitivity  of  the  end-of-tape  markers;  aluminum 


A- 9 


1 


i ’ 


I 


foil  markers  were  placed  to  keep  the  recorder  from  overrunning  the  end 
of  the  tape.  A photocell  sensed  the  reflected  light  from  a small  light 
source,  but  in  field  use  the  intense  natural  light,  even  when  units  were 
shielded,  caused  premature  shut-off. 

20.  The  major  problem  encountered  with  the  recorder  data,  as 
opposed  to  the  actual  recording,  was  with  the  retrieval  of  the  data. 
During  the  recording  of  field  data,  especially  if  12  v-d.c.  was  used  as 
the  power  source,  the  recorders  nay  occasionally  have  skipped  one  or 
more  of  the  16  channels  that  were  being  recorded.  This  presented  few 
difficulties  if  the  location  of  the  lost  channels  were  known.  In  prac- 
tice, this  was  one  of  the  several  drawbacks  of  the  entire  continuous 
recording  system. 

21.  The  format  of  the  data  recorded  on  the  field  units  was  filed 
in  a unique  way  that  was  not  compatible  with  the  normal  data  processing 
computer.  For  this  reason,  it  was  transcribed  from  6.35-mm  data  tape  to 
12.7-mm  computer  tape  in  compatible  format.  This  was  done  by  using  a 
metrodata  tape  transcriber. 


Data  Processing 


22.  After  transcription,  the  field  data  were  processed  to  convert 
the  raw  millivolt  readings  from  the  amplifiers  into  the  absolute  values 
for  each  of  the  parameters. 

23.  The  processing  of  the  recorded  data  proceeded  via  several 
steps.  The  cartridge  tape  records  retrieved  from  the  field  were  refor- 
matted and  transcribed  onto  industrial  standard  seven-track  magnetic 
tapes  using  a Metrodata  model  625  tape  reader  driving  a Digi-Data 
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recorder.  A block  length  of  960  characters  was  used  as  optimal  for 
input  to  the  CDC  6400  computer  used  in  this  study.  This  taped  informa- 
tion was  then  used  as  input  into  a series  of  computer  programs  result- 
ing in  diurnal  curves  of  radiation,  tide,  pH,  and  temperature. 

24.  The  first  program  accessed  was  TRANSCRIBE  (see  Kelly  et  al., 
1975,  for  complete  description),  which  converted  the  recorded  input 
voltages  on  the  seven  track  tape  into  decimal  equivalents.  Simulta- 
neously, it  checked  for  outlier  data  points  and  errors  in  recording. 

The  output  from  this  program  was  written  onto  a permanent  master  tape 
that  could  be  accessed  at  some  future  date. 

25.  At  this  stage  a program,  PROCESS,  was  utilized  to  convert  the 
millivolt  readings  into  absolute  values  (see  Kelly  et  al.,  1975,  for 
detailed  description) . This  conversion  utilized  the  output  tape  pro- 
duced by  TRANSCRIBE  in  addition  to  field  calibration  information. 

Within  the  program,  the  continuously  monitored  information  was  averaged 
for  a period  of  time  specified  by  the  operator  and  then  converted  to 
absolute  values  via  a series  of  mathematical  formulations.  This  infor- 
mation could  be  either  printed,  punched  onto  cards,  or  written  onto 
another  tape  or  disc  file.  For  purposes  of  visual  display  and  analysis, 
the  punched  card  output  was  used  as  input  to  a plot  routine,  which 
plotted  the  radiation,  tide,  temperature,  and  pH  as  a function  of  time. 
This  information  could  be  displayed  in  a variety  of  formats  that  in- 
cluded fixed  and  variable  scales  for  the  y-axis  as  well  as  plotting  all 
environmental  parameters  at  one  station  and/or  comparing  parameters 
between  stations. 
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APPENDIX  C:  SURVEY  SITES 

1.  Plant  and  sediment  surveys  were  conducted  at  a variety  of 
sites  along  the  eastern  seaboard,  from  New  Jersey  to  Georgia,  to 
attempt  to  determine  if  the  pattern  shown  by  the  zonation  of  the  marsh 
plants  was  reflected  by  certain  constituents  within  the  sediments. 

Most  of  these  sites  were  visited  only  once — on  extended  survey  trips — 
although  a subsequent  trip  was  conducted  to  several  of  the  sites  in 
order  to  collect  a biomass  estimate  from  the  end  of  the  growing  season. 
Many  of  the  sites  were  visited  by  students  and  personnel  only  parti- 
ally connected  with  this  study,  and  as  a result,  the  data  sets  are  not 
uniform. 

2.  The  name  and  computer  abbreviation  of  the  survey  sites  are 
listed  in  Table  C-l  and  the  location  of  the  sites  is  shown  in  Figure  C-l. 
The  objective  at  these  sites  was  not  the  creation  of  detailed  vegeta- 
tion maps  and  fine  scale  resolution  of  the  changes  in  chemical  and  bio- 
tic parameters,  but  the  general  depiction  of  trends  in  certain  chemi- 
cal constituents  of  the  sediments,  and  their  relation  to  changes  in 

the  dominant  vegetative  species. 

3.  At  the  survey  sites,  sampling  stations  were  located  on  a 
transect  that  extended  from  the  lowest  elevation  colonized  by  Spartina 
alterniflora.  through  the  S,.  alterniflora  zone  into  the  high  marsh 
where  this  pattern  existed.  The  stations  with  the  lowest  numbers  were 
the  lowest  areas  colonized  by  JS.  alterniflora  and  the  highest  numbers 
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were  generally  the  uppermost  region  of  the  high  marsh,  (Exceptions  to 
this  are  Taskinas  Marsh  and  Bennett's  Marsh,  but  these  are  described 
in  detail.)  Transect  lengths  varied  greatly  but  these,  unfortunately, 
were  not  always  recorded.  Vegetation  samples  and  sediment  cores  were 
attempted  at  each  station.  Single  or  duplicate  samples  of  live  stand- 
ing crop  were  taken;  however,  some  samples  were  lost.  Quadrat  size 
was  varied  according  to  plant  type  and  density;  0.1m2  quadrat  samples 
were  harvested  ir  the  S_.  alterniflora  zone  and  0.04  m2  samples  were 
taken  where  the  much  denser  S^.  patens  or  D.  spicata  predominated. 

These  samples  were  sealed  in  plastic  bags  for  transporting  to  the  lab- 
oratory for  counting  and  weighing.  Table  C-2  gives  the  computer  abbre- 
viations for  the  plant  data. 

4.  Sediment  cores  were  taken  with  either  a 7.6-cm  (3-inch) 
internal  diameter  P.V.C.  corer  or  a 6.35-cm  (3.5-inch  i.d.)  brass  core 
tube.  These  were  extruded  in  the  field,  wrapped  securely  in  plastic 
bags,  and  frozen  with  dry  ice  for  transportation  to  the  laboratory. 

5.  For  certain  samples,  the  pH  and  oxidation-reduction  potential 
(EH)  of  the  sediment  were  measured  immediately  upon  extrusion  of  the 
core.  pH  was  measured  using  a Sensorex  S 200  c combination  pH  elec- 
trode, and  an  Orion  model  401  portable  pH  meter.  Eh  was  measured 
using  the  reference  half  of  the  combination  pH  electrode  and  a Thomas 
4096-D20  platinum  electrode. 

6.  Cores  were  stored  at  -15°C  prior  to  lab  analysis.  After 
thawing  for  12-18  hours  at  room  temperature,  each  core  was  divided 
into  10-cm  sections.  A subsample  weighing  12.5  g (wet  weight)  was 
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cut  from  each  of  these  sections  on  a longitudinal  strip.  These  sub- 


samples, representing  the  three  cores  obtained  from  the  same  sample 

zone  and  depth  range  (a  total  of  37.5  g of  sediment),  were  consolidated 

with  150  ml  of  extractant  solution  (0.025  N HC1  + 0.01  H SO  ) (Nelson 

2 ■* 

et  al. , 1953)  and  homogenized  for  30  seconds  at  high  speed  in  a stain- 
less Waring  blender.  The  resulting  solution  was  filtered  through 
Whatman  #42  filter  paper  and  collected  in  plastic  sample  bottles. 
Several  drops  of  chloroform  were  added  to  each  sample  bottle  to  inhi- 
bit bacterial  activity.  Filtered  samples  were  frozen  until  analysis. 
Cores  from  the  survey  sites  were  treated  similarly,  but  only  one  core 
was  utilized  per  analysis. 

7.  Subsamples  from  all  three  cores  representing  a given  sample 
zone  were  consolidated,  as  above.  Consolidated  subsamples  were  dried 
in  a warming  oven  for  48  hours  at  100°C.  Dried  samples  were  cooled 
and  reweighed,  then  pulverized  with  mortar  and  pestle.  Subsamples  of 
approximately  one  gram  were  transferred  to  crucibles  and  weighed 
within  ± 0.0002  g on  a Sartorius  analytical  balance.  The  samples  were 
then  ashed  in  a muffle  furnace  for  12  hours  at  550°C  and  were 
reweighed  after  cooling  in  a desiccator.  Percent  organic  matter  was 
determined  by  weight  loss  during  ashing. 

8.  In  the  data  section  that  follows,  the  following  units  are 
used:  biomass  (grams  dry  weight/m2),  stem  height  (cm).  All  chemical 
values  are  in  parts  per  million  (p.p.m.)  and  the  day  in  the  Julian 
day  of  the  year. 
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Richmond 


POROPOTANK 


OLLEY  CR. 


CARTER  CR 


MILES 


Location  of  survey  sites  along  the  central  eastern  coast 
of  the  U.  S.  The  two  sites  not  shown  on  the  figure  are 


Avalon,  New  Jersey,  and  Sapelo  Island,  Georgia. 


Table  C-l.  List  of  Sites 


1 

AVA 

Avalon,  New  Jersey 

2 

BAV 

Bavon,  Virginia 

3 

BEN 

Bennett's  Creek,  Virginia 

4 

CAR 

Carters  Creek,  Virginia 

5 

CRA 

Craney  Island,  Virginia 

6 

FIS 

Fishermans  Island,  Virginia 

7 

FOL 

Folley  creek,  Virginia 

8 

GAR 

Gargathy  Inlet,  Virginia 

9 

LYN 

Lynnhaven  Inlet,  Virginia 

10 

POR 

Poropotank  River 

11 

SAP 

Sapelo  Island,  Georgia 

12 

605 

Route  605  Bridge,  Virginia 

13 

TAS 

Taskinas  Creek,  Virginia 

14 

UAL 

Walker  Marsh,  Virginia 
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Table  C-2.  Plant  Identification  Codes 


S.ALT 

D.SPI 

S.PAT 

PHRAG 

SCIR 

J . ROM 

J.GER 

TYPHA 

S.EUR 

S.BIG 

SALIC 

S. DEAD 

T. DEAD 
LITT 
UNK 
SP-DS 
SA-DS 
SA-SP 
P-DS 


Spartina  alternif lora 
Distichlis  spicata 
Spartina  patens 
Phragmites  communis 
Scirpus  sp. 

Juncus  romerianus 
Juncus  gerardi 
Typha  sp. 

Spartina  europea 
Spartina  Bigelowi 
Salicornia  sp. 

Standing  dead 
Total  dead 
Litter 
Unknown 

S.  patens  and  D.  spicata  mixed 
S.  alterniflora  and  D.  spicata 
!S.  alterniflora  and  patens 
£5.  alterniflora,  SL  patens,  and  I).  spicata 
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In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Zieman,  Joseph  C 

Modeling  of  ecological  succession  and  production  in  estuarine 
marshes  / by  Joseph  C.  Zieman,  William  E.  Odum,  University  of 
Virginia,  Department  of  Environmental  Sciences,  Charlottesville, 
Virginia.  Vicksburg,  Miss.  : U.  S.  Waterways  Experiment  Station  ; 
Springfield,  Va.  : available  from  National  Technical  Information 
Service,  1977. 

160,  11,  6 p.  : ill.  ; 27  cm.  (Technical  report  - U.  S.  Army 
Engineer  Waterways  Experiment  Station  j D-77-35) 

Prepared  for  Office,  Chief  of  Engineers,  U.  S.  Army,  Washington, 
D.  C. , under  Contract  No.  DACW39-73-C-0121  (DMRP  Work  L'nit 
No.  4A05) 

Appendix  B on  microfiche  in  pocket. 

References:  p.  155-160. 

1.  Dredged  material.  2.  Ecological  models.  3.  Ecological  suc- 
cession. 4.  Estuarine  ecology.  5.  Plant  growth.  6.  Salt  marshes. 
I.  Odum,  William  E. , joint  author.  II.  United  States.  Army.  Corps 
of  Engineers.  III.  Virginia.  University.  Dept,  of  Environmental 
Sciences.  IV.  Series:  United  States.  Waterways  Experiment  Station, 
Vicksburg,  Miss.  Technical  report  ; D-77-35. 

TA7.W3A  no. D-77-35 
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APPENDIX  B:  MODEL  USER  INFORMATION 


General 


1.  The  model,  MARSH1,  was  written  In  FORTRAN  IV  as  this  Is  the 
most  universally  used  higher  level  scientific  computer  language.  The 
original  design  and  implementation  were  carried  out  on  a Control  Data 
Corporation  6400  machine,  but  the  program  was  written  to  be  machine* 
Independent,  and  machine  or  CDC  software  functions  or  conventions 
which  are  not  available  on  other  brand  computers  were  not  used. 

2.  The  model  structure  was  designed  for  versatility  and  ease  of 
Modification  of  the  input  and  output  data.  Of  particular  note  is  the 
relative  ease  of  using  the  model  from  any  of  a variety  of  interactive 
keyboard  terminals.  If  a computer  syntem  le  available  with  sufficient 
bulk  storage  capacity  to  allow  the  storage  of  the  program,  terminal 
use  is  quite  simple.  The  output  routines  allow  specification  cf  80  or 
132  column  output,  and  other  plotted  output  variables  are  automati- 
cally scaled. 


Data  Input 


3.  The  input  data  necessary  for  the  execution  of  the  model 
falls  into  three  general  categories:  site-specific  da«--.  Includes  the 
numbers  of  stations  along  the  transect,  the  number  and  species  of 
plants  to  be  Included  in  the  parameterization,  and  the  latitude  of  the 
site.  The  model  parameterization  values  include  such  informaclon  as 
the  type  of  growth  algorithm,  the  time  step  for  the  numerical 
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integration,  and  Intervals  for  printing  and  plotting.  The  plant  para- 
meterization section  Includes  plant  specific  intrinsic  growth  rates 
and  the  curves  of  species  response  to  the  driving  variables.  The  most 
frequently  changed  cards,  the  site-specific  data,  are  placed  first  in 
the  data  deck,  in  order  to  facilitate  manipulation. 

4.  The  parameters  associated  with  the  input  data  are  explained 
in  the  next  section  and  should  be  self-explanatory  with  the  exception 
of  those  related  to  the  species  parameterization  curves  and  the  growth 
algorithms.  The  species  parameterization  curves  represent  response  of 
the  plants  to  the  four  environmental  forcing  functions. 

5.  The  x-axls  of  the  parameterization  curve  reflects  the  value 
of  the  specific  forcing  function  (e.g.  temperature),  while  the  y-axls 
represents  a normalized  growth  response  to  this  variable.  The 
response  of  a plant  to  a specific  level  of  an  environmental  forcing 
function  ranges  between  0.0  and  1.0;  if  the  plant  response  is  0.0,  the 
environmental  factor  has  completely  curtailed  metabolic  activity, 
while  a response  of  1.0  Indicates  the  level  of  the  factor  la  optimum 
for  growth. 

6.  The  other  data  input  variable  which  may  require  clarifica- 
tion is  the  one  associated  with  the  growth  algorithm-,  the  mathematical 
method  for  calculating  the  level  of  biomass  for  each  species.  Two 
algorithms  are  available  within  this  program,  the  logistic  and  maximum 
productivity  methods.  Several  assignations  are  inherent  in  using 
either  of  these.  In  the  logistic  formulation  productivity  is  deter- 
mined as  a function  of  the  intrinsic  growth  rate,  the  limiting 


B-2 


B-3 


coefficients,  and  the  contemporaneous  biomass.  Also  included  is  a nega- 
tive feedback  term  associated  with  environmental  Inhibition.  The 
naxiraum  productivity  algorithm  is  the  product  of  the  limiting  coeffi- 
cient multiplied  by  maximum  productivity  per  unit  time.  Either  method 
is  applicable  depending  on  the  assumptions  the  user  desires  to  accept. 

7.  Following  is  a description  of  data  input  and  data  card  pre- 
paration. The  user  should  first  have  read  carefully  Part  VII  of  this 
report  as  this  section  is  designed  to  supplement  that  Part.  The 
example  data  used  here  will  reproduce  the  simulations  shown  in  Figures 
50  to  55. 


Data  Card  Preparation 


I.  General  program  information 


A.  Limitation.)  per  problem 


1.  Number  of  plant  species  < 5 

2.  Number  of  time  intervals  < 500  (months) 


3.  Number  of  sites  in  marsh  < 10 


B.  Output  options  available  from  program 


1.  Graphical  output  of  the  standing  crop  of  each  species,  as 
a function  of  time. 

2.  Graphical  output  of  plant  growth,  as  a function  of  space 
and  time. 

3.  Printed  output  of  the  growth  of  each  species,  as  a function 
of  time. 

II.  Data  Card  Preparation  (each  data  set  contains  one  or  more  data 

cards).  The  columns  to  the  left  of  each  parameter  discussed  indi- 
cate the  location  where  the  Information  should  be  placed  on  the 
data  card.  All  information  typed  on  the  data  cards  should  be 
right-justified.  This  means  that  if  there  are  5 columns  design- 
ated for  a variable  and  it  is  represented  by  2 digits,  place  them 
in  the  4 and  5 columns,  respectively.  In  addition,  a listing  of 
the  data  format  in  FORTRAN  is  included  at  the  beginning  of  each 
data  card  description.  The  actual  numbers  typed  on  the  input 
cards  may  be  either  real  or  integer  (Integer  - no  decimal  point, 
real  has  decimal  point).  This  is  specified  for  each  variable. 


A.  Set  1 (1  card)  FORMAT (4 X.  II.  2X,  13.  IX,  F4.2,  4X,  II). 


col  5 NUM(lnteger)  = the  number  of  species  to  be 

simulated. 

col  8-10  IMAX(integer)  ■ the  maximum  number  of  time 
units  per  simulation  (i.e.  months). 

col  12-15  DT(real)  = the  time  interval  for  the  integra- 
tion routine  (between  0 and  1). 


col  20  METHOD(integer)  ■ variable  designating  the 

productivity  algorithm,  1 « logistic  formu- 
lation, 2 “ maximum  productivity  method. 

col  25  NCASE(integer)  * number  of  cases  of  separate 

simulations. 


Set  2 (3  card  set  for  each  species) 

1.  Card  1 FORMAT (IX,  F4.1,  5X.  F5.0.  IX,  F4.0) 


col  2-5  R(real)  “ Intrinsic  growth  rate  for  the  plant 
under  consideration. 


col  11-15  PEAKB(real)  « peak  standing  crop  for  the 
plant  under  consideration. 

col  17-20  VAL(real)  » initial  biomass  to  start  simula- 
tion for  the  species  under  consideration. 

Card  2 (scaling  factor  (0.0  to  1.0)  uniformly  scales  up 
or  down  response  of  plant  to  an  environmental  parameter) 
FORMAT (4  (2X,  F3.1)) 

col  3-5  LST(real)  ■ scaling  factor  for  light  para- 
meterization curve. 

col  8-10  TSF(real)  » a scaling  factor  for  temperature 
parameterization  curve. 


col  13-15  SSF(real)  = scaling  factor  for  the  salinity 
parameterization  curve. 

col  18-20  HSF(real)  » scaling  factor  for  the  tidal 
Inundation  parameterization  curve. 


3.  Card  3 FORMAT (PA 10) 

col  1-10  ISPEC1E  ■ species  name  (i.e.  S.  alternlf lora). 

Note:  Remember  that  this  sequence  of  three  cards  Is 
repeated  for  each  species. 


C.  Set  3 (general  marsh  data,  1 card)  FORMATS (IX,  F4.0)) 


col 

2-5 

ZLAT  (real)  **  latitude  of  marsh  under  simula- 
tion . 

\ 

r 

col 

7-10 

XSAL  (real)  « average  marsh  creek  salinity  In 
parts  per  thousand. 

t 

col 

12-15 

RANGE  (real)  * mean  tidal  range  in  feet  (source, 
NOAA  tidal  tables). 

i 

col 

16-20 

SPRING  (real)  - spring  tidal  range  In  feet 

(source,  NOAA  tidal  tables). 

D*  Set  4 (general  site  A information,  3 cards) 


"*■*  ' ■■  ■ ••J'V '• 
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If: 

NFLAC  » 0,  plotting  Is  desired. 

IFLAC  » 1,  printing  is  desired. 

IFLAG  « 2,  both  printing  and  plotting 

j is  desired.  , 

f 

I col  10  NFLAC  ■ parameter  specifying  style  of  plotting 

If: 

NFLAC  - 0,  plotting  species  growth 

through  time  at  a particular 
station. 

NFLAC  « 1,  plotting  a species  growth 
through  space  and  time. 

col  14-15  IPLOTDT  - plot  interval  (number  of  tine  inter- 
vals between  plotted  points). 

col  19-20  IRPRINDT  •-*  print  interval  (number  of  time 
t intervals  between  printed  points). 

col  24-25  INTER  “ number  of  blank  spaces  between  plotted 
: data  on  graphic  output. 

2.  Card  2 F0RMAT(4X,I1,5X.  F5.0) 

col  5 KF1AG  ■ parameter  designating  how 

the  y-axis  is  computed. 

If: 

KFLAG  - 0,  all  variables 
using  a combined 

• the  y-axis. 

KFLAG  “ 1,  each  variable 
using  a separate 
y-axis. 

KFLAG  » 2,  all  variables 
using  a combined 
the  y-axis  which 
by  the  operator. 

col  11-15  QMAX  ■ operator  specified  range  for  the  y-axis 
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the  range  on 


are  plotted 
range  for 

is  plotted 
range  on  the 

are  plotted 
range  for 
is  specified 
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Card  2-n.  FORMAT (11( XX. F6. 2)) 

The  three  variables  to  the  right  are  typed  sequentially  on 
as  many  data  cards  as  necessary.  Each  data  point  within 
each  variable  occupies  6 spaces  (right  justified)  and  is 
separated  bv  1 column. 

CLEN  » the  endpoint  of  each  line  segment  (i.e. 

value  on  x-axis  corresponding  to  end  of 
line  segment).  In  the  example  CLEN  » 
400.0,  600.0,  1000. 0. 

LL1  » y-axis  intercept  of  each  line  segment. 

LL1  “ 0.0,  1.0,  2.5  for  each  respective 
line  segment  in  the  above  example.  Note 
that  the  y-intercept  of  segment  3 is  2.5 
even  though  the  parameterization  curve 
never  reaches  this  value. 

LL2  * slope  of  each  line  segment.  In  this 
example  LL2  “ 1.0,  0.0,  -1.0. 

In  this  example  the  2nd  card  would  appear  as  follows  400.  Q 
600.0,  1000.0,  0.0,  1.0  2.5,  1.0,  0.0,  -1.0.  The  paramet- 
erization curves  may  have  any  number  of  segments.  If  the 
curves  require  more  spaces  than  are  available  on  one  card, 
simply  continue  on  with  as  many  cards  as  is  necessary 
using  the  same  format. 

Set  9 (temperature  parameterization) 

The  parameterization  of  sets  9 to  11  (temperature,  salin- 
ity, and  tidal  height)  is  identical  to  the  method  of  para- 
meterization for  solar  radiation. 

Card  1 FORMAT (4X, 12) 

col  5-6  NUML  - the  number  of  line  segments  in  the  para- 
meterization curve. 

Card  2 F0RMAT(11(1X,F6. 2)) 

CTEN  » endpoints  of  each  line  segment. 

TT1  “ y-axis  intercept  of  each  individual  line 
segment . 

TT2  « slope  of  each  line  segment. 


Sample  Data  Set 


8.  The  following  example  data  are  a card-image  reproduction  of  the 
input  cards  needed  to  produce  the  results  seen  in  Figures  50  to  55  of  the 
main  text. 
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Program  Listing 

The  program  listing  is  reproduced  on  pages  B-12  through  B-89. 
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